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Abstract 



In this book, we give a full expository lecture about the Ricci-Hamilton flow on 
surfaces. Hence we derive the Li-Yau-Hamilton's Harnack inequality, the entropy 
estimate of R.Hamilton, the isoperimetric estimate , and the improved gradient 
estimates of W.X.Shi, B.Chow, Bartz-Struwe-Ye, for the Ricci-Hamilton flow. We 
discuss the Maximum principle for tensors of R.Hamilton and the moving plane 
method. We also derive in full detail some estimates for curve shortening flow in a 
Riemannian manifold. We give a detailed computation for our estimates involved. 
We also introduce G.Perelman's W-functional and give some geometric application. 
In an appendix, we consider the Nirenberg problem on S 2 by using the perturbation 
method. 

This book is based on my lectures given at Tsinghua University at Beijing, 
National University of Singapore, and Nankai University at Tianjin given in the past 
three years. My aim for these lectures is to introduce graduate students into this 
exciting area of mathematics. I am sure, with a nice understanding of the material 
presented in this book, people can read works of R.Hamilton and G.Perelman well. 
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Preface 



I made up my mind to write this book in December 2002 just after the appearance 
of G.Perelman's first paper on Ricci-Hamilton flow. It is clear that earlier papers of 
R.Hamilton on Ricci-Hamilton flow are beautiful and key point oriented, however, 
many details and further explanations have to be given for graduate students and 
non-experts to understand. In my lectures in this direction, I often found difficulty 
to explain the detail for the Maximum principle for tensors. So I add some de- 
scription in an appendix. As us know that for the scalar linear partial differential 
equations, the positivity of the first eigenvalue of the corresponding operator implies 
the maximum principle, so one can expect such an result true for tensors. 

Once we have a local existence of the flow, we try to extend the flow. That 
is equivalent to finding the long-term behavior of the system. To study the global 
behavior of the flow, one often likes to do a priori estimates for the flow, then one can 
use the Arzela-Ascoli theorem to get a convergence result. Actually, this principle 
is often used in the study of the Ricci-Hamilton flow. However, the estimates are 
not easy. The guiding line is to use the maximum principle. The first difficulty 
is to choose nice geometric quantities (often scaling-invariant quantities) and find 
nice equations (like the Bochner-type formula). In geometric problems, the Bianchi 
identities are always the conservation laws. Then the symmetry or invariant property 
preserved by the flow helps to find the key point to do analysis for convergence. So 
one has to study the self-similar or soliton problem. The common feature is that the 
soliton equation is a fixed point equation. Again one has to face a priori estimate. 
Occasionally, with the help of symmetry, one can reduce the soliton equation to an 
nice ordinary differential equation and luckily find solitons by solving the dynamical 
system. An important trick in finding the a priori estimate is the blow up analysis 
for the flow. Then one has to meet the limiting flows (sometimes it is a soliton 
problem). One may call such problem as singularity analysis. This is like the tangent 
cone problem in minimal surface theory. Classifying the limits and studying which 
singularities are real ones in the flow are the essential objects. In this program, 
geometric and analytical considerations play important roles. The key part is the 
Harnack type estimate. A beautiful work done by P.Li and S.T.Yau in 1986 serves 
as a model for such goal. R.Hamilton further developed this estimate. By now, such 
an estimate is called Li-Yau-Hamilton Harnack estimate. 

In some cases, one can expect global convergence, so one tries to do pinching 
estimates. Pinching set may be defined by convex or concave functions on matrix 
space. R.Hamilton found that the associated ordinary differential equations ( i.e. 
ODE's for short ) help this goal a lot. Usually, the ODE's are monotone systems. 
They are rare in the dynamical system, but often occur in geometric evolution 
equations. Finding nice flows in geometry with nice pinching set are important. 

All these ideas may be explained in the Ricci-Hamilton flow on closed surfaces. So 
this is a good model to introduce non-experts or graduate students to this fascinating 
mathematical area. This is our aim for writing this book. We have to point out 
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that the geometric computations are easy to follow in this case. The other feature is 
that in dimension two, the Ricci-Hamilton flow and the Yamabe flow are the same, 
so one can easily get the local existence of the flow in this case. 

The structure of this note is as follows. In section two, we prove the local 
existence of the flow. In section three, we derive the evolution equation for curvature 
and prove that the global flow exists. Section four studies the behavior of the solution 
at infinity in the case when the average of the scalar curvature r < 0. In the section 
five, we consider Li-Yau-Hamilton inequality of Harnack type for the curvature. In 
section six, we derive the decay entropy estimate. In sections seven and eight, we 
give the uniform estimates for curvature and its derivatives. In Chapter one, we 
give a full argument of Theorems 1.1.1 and 1.1.2. 

In Chapters two and three, we present two other proofs of Theorem 1.1.2 on the 
sphere based on the works [4] and [19] 

We introduce G.Perelman's entropy functionals in chapter four. We give some 
basic properties of these functionals and geometric applications. 

We add four appendices, and we hope this will help readers' understanding of 
key ideas used in the Ricci-Hamilton flow. 

The lesson we learned from the works of Li-Yau , R.Hamilton and G. Perelman is 
that a deep understanding of geometric parabolic partial differential equations can 
yield important new results in Riemannian Geometry and in Differential Topology. 
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Chapter 1 

Ricci-Hamilton flow on surfaces 



1.1 Introduction 

In 1988, R.Hamilton [17] studied the evolution of a Riemannian metric (g^) on a 
compact surface by its scalar curvature R under the flow 



where r is the average value of R. We will call this flow the normalized Ricci- 
Hamilton flow on M. In short, we just call it the Ricci-Hamilton flow in this chapter 
except expressed explicitly. It is clear that this flow make sense in higher dimensions 
and it is called the Yamabe flow, which has been studied by many people. We will 
call such a flow on surfaces the Ricci-Hamilton flow. 
R.Hamilton [17] proved the following 

Theorem 1.1.1 For any initial data, the solution exists for all time. For the con- 
vergence at infinity, we have that 

(1) If r < 0, then the metric converges to one of constant curvature; 

(2) If R > 0, the metric converges to one of constant curvature 

The global existence of the flow was obtained early by H.D.Cao in [5]. The proof 
of Hamilton's Theorem 1.1 is based on two outstanding estimates. One is the so 
called Li-Yau-Hamilton inequality of Harnack type for the scalar curvature. The 
other is a decay estimate for the entropy which is defined as 



We can use the remarkable work of G. Perelman [30] to study the Ricci Hamilton 
flow on surfaces, but we will not present it in this book. 

Later, B.Chow [9] removed the condition in Theorem 1.1 (2). His result is 

Theorem 1.1.2 For any metric on S 2 , then under the Ricci-Hamilton flow, the 
scalar curvature becomes positive in finite time. 



dt9ij = (r - R)g, 



(1.1.1) 




(1.1.2) 
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Combining the two theorems above yields 

Corollary 1.1.3 // (gij) is any metric on a closed surface, then under the Ricci- 
Hamilton flow, the metric converges to one of constant curvature. 

This last result implies what we called the (differential geometric) uniformazation 
theorem. It says that 

Theorem 1.1.4 any Riemannian metric on a two dimensional closed surface is 
point-wise conformal to a metric of constant curvature. 

The Ricci-Hamilton flow on surfaces is also studied by B.Osgood , R.Phillips and 
P.Sarnack [29] and J.Bartz, M.Struwe, and R.Ye [4], and others. The work of [4] 
used a moving plane method introduced by A. D. Alexandrov [2]. 

We use the following notations for the Riemannian metric surface (M 2 , g) except 
explicitly stated. 

A g is the Laplacian operator of the metric g, and 
D is the covariant derivative of g. 

For a metric family (g(t)), we let 
R be the scalar curvature of g = g(t) 

In following, we write 
Rmax as the maximum of R at time t. 
Rmin as the minimum of R at time t. 
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1.2 preliminary material for Ricci-Hamilton flow 
on surfaces 



Let 

g = gijdx % dx 3 

be the Riemannian metric on the closed surface M, written in local coordinates (x l ). 
Let 

(g ij ) = (9a)- 1 

be the inverse of the matrix (gij). Set 



H = a/ detgij 

be the volume density. Then along the Ricci-Hamilton flow, we compute 

dtn = hr g (g t )fi = ^g lJ d t gijii = (r - R)/i 
Hence for the total area A = f M dfx, we have 

A t = / (r - R)dfi = 0. 
Jm 

That is, the total area is preserved along the flow. 
By the Gauss-Bonnet formula we have 



Rdfx = 4vrx(M) 



M 



where x(M) is the Euler number of the surface M. A well-known fact in Topology 
is that 

X (M) = 2(l-g), 

where g is the genus of the closed surface M and it is the number of "holes" of M. 
For example, if M = T 2 is the torus, we have g — 1, and then x(M) = 0. We recall 
that the usual way to write the Gauss-Bonnet formula is the following 



/ Kdii = 27r X (M) 
Jm 



where K — R/2 is the Gauss curvature of M. 
Since 



/ Rdfi/A, 
Jm 



we get that 

r = 4irx(M)/A 

which is a constant along the flow, furthermore, if x(M) is negative, we have r < 0. 
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1.3 Short-time existence of the flow 

We point out that the flow exists at least in short time. We will prove this fact by 
using the conformal transformation. In fact, write g = e u g , then we have 

R = e- u (-Au + R ) 

(see appendix for derivation for this formula) and the Ricci-Hamilton flow is reduced 
into the following non-linear parabolic equation 

ut = e- u (Au- R + re u ), (1.3.1) 

with initial data 

u\t=o = 1- 

Many people just call this equation the Yamabe flow. 

By the standard linearization and fixed point argument we can always find a 
short-time solution of (1.3.1) for any initial metric go- Hence we have 

Theorem 1.3.1 For any metric on a compact 2-manifold M 2 , then there is a pos- 
itive constant T such that the Ricci-Hamilton flow exists at the time interval [0, T). 

We now consider the energy estimate for this Yamabe type flow. 
Define 

E{u)= [ (hD u\ 2 + R u)alfi 
Jm * 

Then its first variation is 

5E{u) = {R- r)g, 
Hence along the Ricci-Haimiton flow, we have 

E t = - I (R- rfdfxg. 
Jm 

Clearly from this we have the following estimate 

I alt [ (R g -r) 2 dn < E(u ). 
Jo Jm 

where T is the maximal time for the flow to exist. The second variation of E(-) is 

5 2 (u)(v)= [ \D v\ 2 dfx >0 
Jm 

Hence we can expect the global existence of the flow and its convergence at infinity. 

The following Harnack type estimate is a key for the flow to convergence at 
infinity. 
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Assertion 1.3.2 Assume M = S 2 . Then along the Yamabe type flow, we have a 
uniform constant C > such that 

inf e u > Csupe". 

* t 

We will call this assertion as Bartz-Struwe-Ye estimate since they obtained this 
estimate in [4]. This assertion will be proved in later sections. If we have this 
assertion, then we can use the volume restriction 




to get the uniform L°° bound on u. Using the standard linear theory of parabolic 
equations, we can get uniform bounds for the higher order derivatives of u. Then 
we can extend the solution in short time to a global one. In fact, if T* < +oo is the 
maximal existence time for the Yamabe type flow, we shoud have 

lim u max = +oo. 

This is a contradiction to our uniform L°° bound on u. Here, we point out that all 
metrics g(t) = e u ^go are uniformly equivalent to the metric g . Recall that our flow 
can also be written as u t = r — R. From this, we can get the uniform bound of all 
derivatives of R. 

Using the energy bound 

/ dt \ Ut \ 2 d[x= / dt (R g -r) 2 dfx< E(u )., 
Jo Jm Jo J m 

we can find a sub-sequence tj — > +oo such that u t (tj) converges to zero and Rj = 
R{tj) converges to r in L 2 (M). Using the Arzela-Ascoli theorem, we have that u t (tj) 
converges to zero and Rj = R(tj) converges to r in C k (M) for any k > 2. In later 
sections, we can see that the whole g(t) converges to a constant curvature metric. 

We also point out that the Ricci-Hamilton flow is not strictly parabolic system, 
so one can not conclude the short-time existence of the flow from the standard 
argument. In fact, R.Hamilton [18] first used Nash-Moser implicit function theorem 
to obtain such a existence. Later De Turck [13] gave a short and beautiful proof 
based on a normalization argument. Their arguments work for any dimensions. 
However, in dimension two, we can avoid all these difficulties by using a conformal 
deformation trick as above. 
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1.4 Global existence of the flow 

We now study long-time behavior of the normalized Ricci flow on the compact 
surface M: 

dt9ij = (r - R)gij. 

Then using the formula 

R = e- u (-Au + Rq) 
we have the flow for scalar curvature as follows. 

d t R = AR + R(R-r). (1.4.1) 

Applying the maximum principle to the equation (1.4.1) we know that if the scalar 
curvature at t = is non-positive, it remains so for t > 0, and if the scalar curvature 
at t = is non- negative, it remains so for t > 0. Clearly, we have the following 

Proposition 1.4.1 There is a uniform constant C such that along the flow, we 
have 

R> -C. 

In fact, we have the following two assertions: 

(1) . If r > and the minimum of R at initial time is negative, the minimum of 
R is increases. 

(2) . If r < and the minimum of R at initial time is less than r, the minimum 
of R is also increases. 
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1.4.1 negative scalar curvature 

Using the maximum principle we have the following result 

Theorem 1.4.2 Assume that there are two positive constants C2 > C\ > such 
that —C2 < R < —C\ at t = 0. Then the scalar curvature R remains so and there 
is a constant C > such that 

re -c 2 t < r _ R < Ce rt 

So the flow exists for all t > and R approached r exponentially as t — > +00. 

Proof: Note first that by the definition of r, we have — C 2 < r < —C\. Then for 
fixed t > 0, at the maximum point of R we have R — r > and the differential 
inequality 

^R < R(R-r) < 0. 
at 

Hence we have R < —C\. Using this and the differential inequality again, we have 

^R < R(R - r) < -CAR - r). 
at 

So 

(log(R-r))'< -d. 
and by integrating, we find that 

i?-r < (R(0) - r)e- Clt 

and 

r-R>(r- R(0))e~ Clt > re^ . 
At the minimum point of R, we have 

^-R > R(R-r) > r(R-r). 
Then we immediately obtain that 

r-R< C 2 e rt . 
Therefore, we have proven the theorem . □ 
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1.4.2 positive scalar curvature 

For the positive scalar curvature case, the flow is hard to treat. In this case, we have 
the following 

Theorem 1.4.3 . Assume that r > and R is not a constant at t = . 

(1). If R/r > c> at t = 0, then c < 1 and along the flow, we have the estimate 

„ rc 
R > 



c+ (1 - c)e rt ' 



(2). If there is a positive constant C such that R/r < C at t = 0, then we have 
that C > 1, and for t < £ log we have 



R < Cr 



rt 



C-(C- l)e 

Proof: . Since R(0) > 0, by the maximum principle, we have R > along the flow. 

In the case (1), we have at the minimum point of R, denote by / := R m in the 
minimum of R, we have 

1 > Rmin/r > C 

Since at time t — 0, R is not a constant, we have R m in < r, and we get c < 1. Note 
that for t > 0, at the minimum point of R, we have the differential inequality 

l; t ! >Hf-r). 

Hence we have 

-(r^i-r/- 1 . 

Integrating this we find 

f-R > TC 



c+ (1 - c)e rt ' 

By the same method we can prove that C > 1 in the case (2) and the estimate 

R< Cr 

~ C-(C-l)e rt ' 

Then we are done. □ 
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1.4.3 Scalar curvature positive somewhere 

In this part, we assume that the scalar curvature R at t = is positive somewhere 
on M. In this case, we introduce a potential function /. 

Definition 1.4.4 Let (M,g) be a compact surface with a riemannian metric g. The 
potential function f is the solution of the Poisson equation 

Af = R-r 

with mean value zero, that is, f M fd/i = 0. 

We have to consider the existence of such a potential. In fact, we can introduce 
the variational functional 

J(u) = [ (\Du\ 2 /2 - (R-r)u) 
Jm 

on the space 

A:={H\M); [ fdfi = 0}, 

where D is the covariant derivative of g. Using Poincare's inequality we can easily 
see that 

inf J(u) > — oo. 
A 

Then one can minimize J(-) over the Hilbert space A and shows that the minimum 
is attained at some point in A. Obviously, the minimum point of ■/(•) is a solution 
of the Poisson equation 

Af — R — r. 

By the maximum principle, we know that the solution with zero mean value property 
is unique. For more detail, see Th.Aubin's book [3]. 

We now compute the evolution equation for the potential function /. 

Proposition 1.4.5 Let g = g(t) where g(t) is a solution of the Ricci-Hamilton flow. 
Then the potential function satisfies the following evolution equation: 

f t = Af + rf-b 

on M, where 




is only a function of the time variable t. 
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Proof: We do computation. 
Note that 



A(A/ + r/) = A(R-r + rf) 
= AR + rAf 
= AR + r(R-r) 



= Rt + rR- R 2 + rR-r 2 
= R t -(R-r) 2 

We also have 

Rt = d t Af 

Since in local coordinate, we have 
We will write 

f . d 'f 

for short notation in this part. Notice that g^dfT^ = 0. Then we have 

d t Af = g yf l3 + Af t = Af + {R-r) 2 . 

By this we obtain that 

A(A/ + r/) = A/ t . 

This implies that 

Af + rf = Af t + b 
for some constant function b only in space variables. In other word, we have 

Af t = Af + rf-b. 

Differentiating the relation 

/ fdfi = 0, 

J M 

we get 

/ (f t + f(r-R))diM = 0. 



By this we conclude that 



f (rf-b-fA)d» = 0, 
Jm 



which clearly implies that 

r \Df\ 2 



□ 



M 
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Definition 1.4.6 We define a P-function 



h = Af+\Df\ 2 

and a trace-free part of the second order derivative of the potential f: 

Mij = DiDjf - ^A/ • gij. 

We remark that by definition, a P-function is a smooth function which can be 
applied by the maximum principle. As pointed out by R.Hamilton [17], once we 
have a nice bound for h, we can bound R because of the formula 

R = h- \DF\ 2 + r. 

Note also that 

h = e- f div(e f Df). 

Proposition 1.4.7 Along the Ricci- Hamilton flow, we have 
(!)• 

ht = Ah-2\M\ 2 + rh 

where at the normal coordinate at a given point, \M\ 2 = M 2 j. 
(2). 

(d t - A)\M\ 2 = -2\DM\ 2 - 2R\M\ 2 . 

Here we just compute the equation for h and we postpone the proof for (2) in 
the later section (see Lemma 1.9.2). 

Proof: We do computation. First, we have 

\M\ 2 = \M l0 \ 2 = \D l D ] f\ 2 - l -{R-r) 2 . 

Secondly we have 

d t h = d t Af + 2Df-Df t + (R-r)\Df\2 

= Af t + (R-r)2 + 2Df-Df t + (R-r)\Df\2. 

Thirdly we have 

Ah = A(Af+\Df\ 2 ) 

= A(Af) + A\Df\ 2 

= A(Af) + 2\D 2 f\ 2 + 2ADf-Df 

= A(Af) + 2\D 2 f\ 2 + 2DAf-Df + R\Df\ 2 . 
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Finally we have 



(d t -A)h = A{f t -Af) + {R-r) 2 -2\D 2 f\ 2 
+ 2Df-D{f t -Af)-r\Df\ 2 

= A (r/ -b) + (R- rf - 2\D 2 f\ 2 + 2Df ■ D(rf - b)) - r\Df\ 2 
= rAf + (R-r) 2 -2\D 2 f\ 2 + r\Df\ 2 
= r(R-r)(R-r) 2 - 2\D 2 f\ 2 + r\Df\ 2 
= (R-r) 2 -2\D 2 f\ 2 + rh 
= -2\M\ 2 + rh. 

□ 

Applying the maximum principle to the equation 

ht = Ah-2\M\ 2 + rh 

we immediately find that 

dth max ^ T'h max , 

where h max is the maximum value at time t. And this implies that there is a constant 
C > such that 

h < Ce rt . 

Importantly we have 

R = h- \Df\ 2 + r<h + r< Ce rt + r. 

This gives us a upper bound on R for all t. Combining this with Theorem 1.4.3 we 
have the following estimate 



Theorem 1.4.8 For any initial riemannian metric on the compact surface M , there 
is a uniform constant C such that along the flow, we have 

-C <R< Ce rt + r. 

Furthermore, if r < 0, then the scalar curvature remains bounded both above and 
below. 



Using this Theorem and the short-time existence of the flow, we can obtain a 
global solution of the Ricci-Hamilton flow. This follows from the standard continuity 
method for heat flow. So we conclude 

Theorem 1.4.9 For any initial riemannian metric on the compact surface M , the 
Ricci-Hamilton flow has a solution for all time. Furthermore, if r < 0, the global 
solution converges at infinity to a metric with constant negative curvature. 
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1.4.4 average scalar curvature flat case 

Basically we just give some remarks in this part. 

Assume that r = at t = 0. We solve the Poisson equation 

A u = Rq, on M. 

Here Ao and Ro are the Laplacian operator and scalar curvature of the initial metric 
go. It is always solvable since r(0) = 0, where 

r(0) = / R dfi/A. 

Set 

9 = e u g . 

Then we have 

R = e- u (-A g() +R )) = 0. 

In this case, g is a stable point of the Ricci-Hamilton flow. So we may expect that 
the Ricci-Hamilton flow converges to a flat metric at infinity. Using this g to replace 
the initial metric go, we restrict the metric family (g(t) in the conformal class of the 
initial metric g to solve the Ricci-Hamilton flow. Write 

git) = e u g. 

Then the scalar curvature of g(t) is 

R = e- u (-A s u) 
and the Ricci-Hamilton flow is reduced into 

u t = e- u (A- g u) (1.4.2) 
In the computation of this section, we will always write 

A = A~ g D — Dg. 

and dp, = d\i g . 

Using the maximum principle we get that there is a uniform constant C > 
such that \u\ < C . Therefore, the metrics g(t) are uniformly equivalent for all time 
t. This gives us the uniform control of diameter, injectivity radius, and the best 
constant in the Sobolev inequality. 

Theorem 1.4.10 The Ricci-Hamilton flow exists for all time, and it converges at 
infinity to a flat metric. 
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Proof: Using the Ricci flow we have 

J \Du\ 2 dfl + 2 J e- u (Au) 2 dfi = 0, (*) 



d 
di 



and using the bounded-ness of u and the inequality 

J (Aufdft >cj \Du\ 2 dfi 
(a proof of this inequality can be obtained by arguing by contradiction), we have 

j t J \Du\ 2 dfi + c J \Du\ 2 dfi<0 
By the Gronwall inequality we have 

\Du\ 2 djj < Ce~ ct . 
Going back to (*) again and integrating it over time we find 
2 J dt J e- u (Au) 2 dJi < J \Du\ 2 dft(T) 

which implies that 

r+oo p 

2 J dt J R 2 dfi < Ce~ cT . 
By this we can have at least one point t e [T, T + 1] such that 

R 2 dfx < Ce~ cT . 



d 
dl 



On the other hand, using the evolution of R we have 

R 2 dfi + 2 J \DR\ 2 dfi = J R 3 dfx. (**) 



Since R is uniformly bounded, we have 

R 2 dfi < C I R 2 dfi. 
Using the Gronwall inequality again we have 



d 
dt 



R 2 dfi(t) <C J R 2 dfi(i), 
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for any t G [T, T + 1] , which clearly implies that 

R 2 djj < Ce- Ct . 



all t. That is to say that R — > at infinity in L P (M, g ) for any p > 1. 
Integrating (**) we have 



rft y |£>i2| 2 dji < e~ cT , 
and arguing as before we find some t in each [T, T + 1] such that 

J \DR\ 2 dfi < Ce~ ct . 



Note that 

d 



dt 

Using the Holder inequality we have 



DR\ 2 dfi + 2 / \AR\ 2 dfi < -2 / R 2 ARdfi. 



2| J R 2 ARdfx\ < J (AR) 2 dfi + J R 4 dfi. 



d J \DR\ 2 dfi + J \AR\ 2 dfi< J R 4 dfi. 



By this we have 

dt 

Using the inequality 

J \AR\ 2 dfi >cj \DR\ 2 dfx 

we have 

j t J \DR\ 2 dfx + c J \DR\ 2 dfi< J R A dfx. 
Using this we can get 

\DR\ 2 dfi < Ce- Ct 



for some uniform constant c and for all t. Again integrating over time we have 

dt [ (AR) 2 dfx < Ce~ cT 



so we have 

J (AR) 2 dfi < Ce~ ct 

at least in some point in each [T, T + 1] . Compute the evolution of AR we can 
the bound 

,-ct 



J (AR) 2 d/i < Ce 
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for all t. Using this bound, we can bound the L 2 norm of D 2 R. In fact, we have 

J \D 2 R\ 2 dfx = J (AR) 2 -^J R\DR\ 2 dfx. 
Using the L°° bound of u, we have 

\D 2 R\ 2 dii < Ce~ ct 



I 



We now can bound the maximum of R by the L 2 norm of R, DR, and D 2 R. 
Thus we get that R converges to zero exponentially. 

Using the formula u t = —R and the L°° bound of u again, we know that u 
converges to a positive constant at infinity. □ 

To conclude this section we make a remark. To understand the equation (4.1) 
better, we need to study the corresponding ODE: 

s t = s(s - r). 

It is easy to see that the unique solution is of the form: 

where c> 1. Note that, for r > 0, s(t) < 0. 

From the equations for R and s we conclude that 

3 t (R -s) = A(R -s) + (R-s)(R-r + s). (1.4.3) 

By the maximum principle, it follows that if s(0) < minMR(x,0), then R — s > 
for all time. 
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1.5 Li-Yau-Hamilton's Harnack estimates 



1.5.1 case one:R min > 

Set L = logi? and define 

Q = AL + R-r 

We will call Q the Harnack quantity for our Ricci-Hamilton flow. Then we compute 

Rt — RLt, 



and 



Lt ~ R 



Hence, we have 

Using 
we compute 



AR + R(R - r) 
R 

AL + |VL| 2 + R — r 



Q = L t -\VL\ 2 
Q = AL + R-r 



Q t = (R — r)AL + AL t + RL t 

= AQ + A\VL\ 2 

+ (R-r)AL + R(AL + \VL\ 2 + R-r) 

= AQ + 2((VAL,VL) + D 2 L\ 2 ) 

+ 2R\VL\ 2 + (2R - r)AL + R(R - r) 

= AQ + 2(VQ,VL) +2\D 2 L\ 2 

+ 2(R - r)AL + (R- r) 2 + rQ 

= AQ + 2(VQ, VL) + 2\D 2 L + \{R- r)g\ 2 + rQ. 
Therefore, we obtain the following 

Lemma 1.5.1 Under the Ricci-Hamilton flow with R > 0, we have 

Qt>AQ + 2(VQ, VL) + Q 2 + rQ. 
Proof: Using the Hadamard inequality we have 

2|T| 2 > itr g T) 2 

for any symmetric 2-tensor on the surface M. Then we get the inequality from the 
equation above. □ 
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Consider the ODE 



St = s 2 + rs 



with initial data s = s < — r < at t = 0. Solving this ODE we get 

Cr 

s(t) = 

where C = s /(s + r) > 1. Applying the maximum principle we get 

Q(t) > s(t) 

for all t > 0. Hence, we obtain the following Li-Yau-Hamilton Harnack estimate: 

d t \ogR- |Vlogi?| 2 > s(t). 

This estimate is also called the differential Harnack inequality. 

Take two points (xi,ti) and (x 2 , t 2 ) in the space-time, where t\ < t 2 , Xi, x 2 G M. 

Let 

7 : [h,t 2 ] -> M 2 

be a smooth curve joining point X\ and rr 2 . Define 

f-*2 



rt2 

A{x\ti, x 2 t 2 ) — inf / | 7 '| 2 <it 
7 7ti 



Then we have 



2 (|logi?( 7 (0,i) + (Vlogit:,7 / »^ 

"to 



> / 2 (|Vlogi?| 2 - s(f0 - |Vlo gj R| • \i\)dt 
Jti 

> f\sit) -\{\lf)dt 



1 f t2 , /l2 , , Ce rt2 - 1 
|Y | 2 <it — log ■ 



4 / ' " ° Ce rf i - 1 



Therefore, we get 



, 1 x , Ce rt2 - 1 

lo S^7 77 - -7^(^1^2*2) - log -r-r -. 

R(x 2 , t 2 ) 4 Ce rtl — 1 

It is easy to see that there is a uniform constant B > such that 

Ce rt ' 2 — 1 

So we conclude 
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Theorem 1.5.2 Along the Ricci-Hamilton flow with R > 0, there is a uniform 
constant C > such that for all x±,X2 € M and t\ < t 2 we have 

R{ x l,tl) > e -\A{x 1 t 1 ,X2t 2 )-C(t 2 -t 1 ) 

R(x 2 ,t 2 ) ~ 
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1.5.2 case two:R min < 

Following the method of B.Chow [9], we only assume that R > only somewhere. 
We make the choice s = when R > on all M. 

Set L = log(i? — s). We have the following important formula: 

L t = AL + \VL\ 2 + R-r + s. (5.1) 

We now define the Harnack quantity 

Q = L t - |VL| 2 - s = AL + R - r. 

Compute 

d t Q = AL t + (A) t L + R t 

= A(AL + \DL\ 2 + R) + (R- r)AL + AR + R(R - r) 

+ AQ + 2\D 2 L\ 2 + 2(DL,DAL) 

+ R\DL\ 2 + (R — r)AL 

+ AL + \DL\ 2 + R(R — r) 

= AQ + 2(DL, DQ) + 2\D 2 L\ 2 

+ 2(R - r) AL + (R - rf + (R - rf 

+ s\DL\ 2 + (r - s)AL + r(R-r). 

Hence, we have 

Q t = AQ + 2(DL,DQ) +2\D 2 L + ^{R - r)g\ 2 + (r - s)Q 
+ s\DL\ 2 + s(R-r). 

Using the upper bound R < e rt , we have sR > —C. Thus using the inequality 

\D 2 L + 1 -{R- r)g\ 2 > \{AL + R - r) 2 = l -Q 2 , 
we have the following inequality 

Qt > AQ + 2(DL, DQ) + 2\D 2 L +\{R- r)g\ 2 + (r - s)Q + s\DL\ 2 - C 
Note that 

(sL) t = A(sL) + s\DL\ 2 + s(R-r + s) + s(s - r)L. 

Since we have 

L > —C — Ct, 

we obtain that 

(sL) t > A(sL) + 2(DL, D(sL)} - s\DL\ 2 - C. 
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Define 

P = Q + sL 

Then we have 

P t > AP + 2(DL, DP) + Q 2 + (r - s)Q - C. 

Since sL = P — Q is bounded, we have foe some constant C > such that for i 
large enough, we have 

P t > AP + 2(DL, DP) + -(P 2 - C 2 ). 
Applying the maximum principle we have 

1 — ce 0i 

for some c > 1. Therefore, for t large enough, we have 

P > -2C 

and 

Q > SC. 

This is our Li-Yau-Hamilton Harnack estimate. 

We now give an application of this Harnack estimate. Define, for two point (xiti) 
and (2:2*2) in our space-time, 

t 2 



h(x 1 t 1 ,x 2 t 2 ) =inf / \l\t)\ 2 dt 
1 Ju 



where 7 is any path in M connecting X\ and x 2 . Note that, for t 2 > t±, 
L(xiti) - L{x 2 t 2 ) = jf" ^L(-y(t),t)dt = jf" (L t + (PL, 7 '))<ft 

Since 

L t - \DL\2 > -C{t 2 - ti) - L/4, 

Integrating we get 

P(x 2 t 2 ) - s(t 2 ) > e h/4 - c(t2 - h \R(x 1 t 1 ) - 3(h)). 
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1.6 Hamilton's Entropy estimate 

In this section, we assume R > somewhere and we study the following important 
quantity 

R log Rdfi 

which is called the entropy for the Ricci-Hamilton flow, and we will show that this 
function is decreasing in t. 



A 
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1.6.1 case one:R min > 

Following R.Hamilton [17], we define 

Z = r- 1 / QRdfi, 
Jm 

which is a global function in time variable t. Then by the definition of Q, it is easy 
to see that 

Z = r- 1 J(- l ^jf + R(R-r)). 

Compute 

2 

+ R(R - r))(r - R)dn 



i2 

Using the equation 

R t = AR + R 2 - rR 

we find that 

Z t > Z 2 + rZ. 

If Z is positive at some time t , then Z is blow up after a finite time. This is 
impossible by our definition of Z, which always exists along the flow. Then we have 
that Z < 0. Using the expression of Z we find that 

12 



f(-^§- + W-r))<0. 



This means that 



\WR 
~ R 



2 

- > / {R-r f 



Note that 



jR\ogRdfi = J(R-r) 2 -J 



dt J r J y J R 

Then we have 

Theorem 1.6.1 Assume R > at initial time t = 0. T/ien a/ong i/ie global Ricci- 
Hamilton flow, we have that 

R log i?<i/i 



zs decreasing. In particular, we have a uniform constant C > such that 

RlogRdfi < C. 
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1.6.2 case two:R min < 

Note that 

s t = s(s - r). 

So, for v — l/(r — s), we have 

v t = — vs 

We compute and get 

(Rdfi) t = ARdfi 

and 

(sdfj) t = s(s — r)d/i. 

Hence, we have 

[v(R - s)dfj] t = v(AR + s(r - R))dfi. 

This implies that 



d t (v J (R- s) \og(R - s)dfij 

-v [ d/j, + v [ (R — s)(R — r + s)d/i 

J M R — S J M 



< v (R-r) 2 dfi 

JM 

Integrating this inequality in time we find 

p poo p 

/ (R-s) hg(R -s)<C dt {R-rf + C. 

JM JO JM 

Hence, using the estimate in section 2, we get 

Lemma 1.6.2 Along the flow we have a uniform constant C > such that 

[ (R-s) \og(R -s)<C. 

JM 
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1.7 Boundedness of R 

Following R.Hamilton [17], we first study the case when R > in all M 



1.7.1 case one:R > in all M 

Take a point (x,t), where r > 1, in the space-time such that the curvature R at 
(x, r) is the maximum of R at time r. Let 

T = r + R max (r)- l /2 

Then at time T, using the evolution equation of R we have 

dtRmax — Rmax 

Then we have 

— (V-Rmax)t < 1- 

This implies that for t £ [r, T] , 

R m \ x {r) - R m \x(t) <T-r = J R max (r)- 1 /2. 

This gives us that 

Rmaxiy) — ^Rmax{ T ) ■ 

Note, on [r, T], since 

dtftj = (r - < rgtj 

and 

max ( T )9ij 

we know that for any tangential vector X, g(X, X) will grow at most by a constant 
factor. 

Let dT(x,y) be the geodesic distance between the points x and y at time T. 
Then we have 

L(xT,yT) <Cd T (x,y) 2 /(T-r). 
Then by the Li-Yau-Hamilton's Harnack inequality we have 

R(x,r)<R(y,T) 

for all y G B p (x) where p = ix ^R m ax{T)/^- 

By theorem 5.9 in [7] we know that the injectivity radius of M is at least p at 
time T. In the ball B of radius p at x we have that the scalar curvature R is always 
comparable to R ma x(j) > \Rmax(T). Hence we have the estimate 

/ R log Rdp > cp 2 R max (T) log R max (T) = c log R max (T) 
Jb 
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at time T for some uniform constant c > 0. Since 

zlogz + e _1 > 

for all z G R, we have 




Then the entropy estimate show that R max (T) is bounded, and then R m ax(j) is 
bounded. This implies that for all r > 1, R is bounded. 

Using this bound on R, we can bound from below the injectivity radius. Since 
the volume is bounded, this also gives us an upper bound on diameter. Using 
the diameter bound and the fact that the growth of distances is bounded, we can 
conclude that for T — t < 1 the estimate 

L(xr,yT)<C/(T-r). 

By the Harnack inequality again we get for all t > 1, any two point x and y, 

R(x,t) < CR{y,t + l). 

Therefore we get a lower uniform bound of the scalar function R. 
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1.8 derivative estimates of R 

Once we have a uniform C° bound of R, we can get the uniform bounds for all 
derivatives of R for all time. In fact, since we have a uniform bound on volume, 
diameter and the injectivity radius, we have the uniform control of the Sobolev 
constant. So one may use the energy bound method to proceed as follows. 

By induction, we have the following evolution equation for the derivative \D n R\ 2 
of R: 

(d t - A)\D n R\ 2 = -2\D n+1 R\ 2 + E i+j=n D i R*D j R*D n R. 
Integrating over M, we can find some uniform constant n such that 



d 
dt 



\D n R\ 2 + j 2\D n+l R\ 2 < C n J \D n R\ 2 . 
Using the well-known interpolation inequality we have that 



J \D n R\ 2 < C^ 1 J \D n+1 R\ 2 + C J R 2 . 



Then we have 

d 



D n R\ 2 + / \D n R\ 2 < C. 



dt 

By the Gronwall inequality we have that 



I 



D n R\ 2 < C(n), 



for all n and all time. By Sobolev imbedding theorem, we then have the uniform 
bound of all derivatives of R. 

In the interesting paper, W.X.Shi [32] find the following Bernstein type estimate 

Theorem 1.8.1 Along the Ricci- Hamilton flow, if R is uniformly bound, so are the 
derivatives of R. In fact, If \R\ < C , we have, on < t < C~ x , for each > 1, there 
exists a uniform constant C n > such that 

\D n R\ 2 < C n C/t n/2 . 

Proof: Recall the following equations for the scalar curvature and its derivatives: 

(d t - A)R = R(R-r). 

and 

(d t - A)\D n R\ 2 = -2\D n+1 R\ 2 + Y, i+j=n D i R*D j R*D n R. 
Assume that \R\ < C and consider tC < 1. Define 

F = t\DR\ 2 + AR 2 . 
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Then we have 

(d t - A)F < (Ct\R\ - 2A)\DR\ 2 + C 2 A 
Take A > C. We have 

(dt - A)F < C 2 A 

Since F < AC 2 at the initial time t = 0. By the maximum principle, we have 

F < (1 + Ct)AC 2 < 2AC 2 

This gives us 

t\DR\ 2 < 2 AC 2 . 

By induction we can get that there exists a constant C n > such that 

\D n R\ 2 < C n C/t n/2 . 

□ 
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1.8.1 case two:R > only in some where 

In this case we first try to get a uniform upper bound of the scalar curvature function 
R. According to the lesson learned from the case R > 0, we need only to need to 
bound the injectivity radius of g{t). In fact, once we have injectivity radius bound, 
we can combine the Harnack inequality with entropy estimate to get a uniform 
supremum bound of R. The injectivity radius bound also implies a uniform bound 
for the diameter of g(t). From the diameter bound we can get the uniform positive 
lower bound of R — s. Then we can see that R > at large time. Then we obtain 
the following result of B.Chow [9] 

Theorem 1.8.2 Assume that R(0) positive somewhere. Along the Ricci- Hamilton 
flow, we have for large t > 0, R(t) > 0. 



So we need only to prove the following 



Theorem 1.8.3 Along the Ricci- Hamilton flow, we have for any t > 0, we have the 
following injectivity radius bound: 

i(g(t)) > min{%(0)), min J }, 

where i(g(t)) is the injectivity radius of g{t). 

To prove this result, we need the following three lemmas. 

Lemma 1.8.4 Let 7 be a shortest closed geodesic on (M,g). If the length 0/7 is 
less than t= , then 7 is stable. That is to say, the second variation of the length 

V -ft-max 

functional at 7 is non-negative. 

Proof: We argue by contradiction. Assume 7 is unstable. Choose a closed curve 7 
near to 7 such that 

J(7) < Ki). 

By choosing two points p and q in 7, we break 7 into two piece 71 and 7 2 of equal 
length. We can also choose two points p and q in 7, and break 7 into two piece 71 
and 72 of with their lengths < j= We can find unique geodesies fa ( which is 

V ■"■ max 

near by 7$) joining p and q. Then we have 

i(A) < Kid- 

Then we should have 

Assertion: There exists a smooth closed geodesic (3 with 

KP) < Kfaufa). 
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In fact, let A be the space of all non-degenerate pairs (aia^) joining two points 
x G M and y G M with 

l{oti) < 

for i — 1,2. Then A is an open 4-manifold locally parametrized by points in M x M. 
Define 

m = inf(/(«i) + /(a^))- 

Clearly, by using the fact that the exponential map of M being a diffeomorphism, 
we have 

m > i(g). 

As in [28], we can find that the infimum m is achieved by some (5 G A. If /3 is not 
smooth, we would shorten (3 inside A. Hence, we have that (3 is a smooth closed 
geodesic curve. This proves the Assertion. Using this Assertion, we get a contraction 
with our assumption at the beginning of the proof. 
□ 

Lemma 1.8.5 Let 7 be a geodesic loop in g(t) as in the above Lemma. Then along 
the Ricci- Hamilton flow, we have, at time t, 

j t Kl) > rZ( 7 ) 

Proof: . Let T be the unit tangent vector fiend on 7 and let N be its unit normal. 
Note 

V T iV = 0. 

Then we have the second variational for the length functional at 7: 

j^(R{N,T)N,T) = - 1 - Jr>0. 

Then we have 



Hence 

= ! /(r-*)> £(7). 

□ 

Lemma 1.8.6 Suppose % is a shortest closed geodesic in the metric g(t). Ifl{ r i) < 
-, 27T , then for e > small enough, there exists a geodesic 7t _ e in (g(t — e)) with 

Y K ma x(T) 

Klt-e) < Z(7), 

where li^t-e) is the length of 'jt-e ^ n the metric (g(t — e)). 
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Proof: For notation more clear, we write 7 = 7*. Using the lemma above we have 
some e > small such that 

k-e{lt) < kilt), 

where l t - t is the length functional in the metric (g(t — e)). By choosing two points p 
and q, we break r y t into two piece 71 and 72 of equal length in the metric (g(t — e)). 
Since 

k-e(li) < 1 

V K max {t — e) 

for sufficient small e, we can find unique geodesies fa ( which is near by 7$) joining 
p and q in the metric (g(t — e)). Then we have 

Z t _ e (/3i) < Z t _ e (7i) 

Using the Assertion above we have a smooth closed geodesic /3 such that 

Z t _ e (/3) < k-e(faUfa) < k(lt), 

Then we are done. □ 

Proof: According to the lemma above, along the Ricci-Hamilton flow, if the length 
functional of shortest closed geodesic is less than , 2n , it is increasing. Then 

V Kmax (t ) 

by using the Klingerberg lemma [7] , we have 



7T 

i(^(t)) > min{i(g(0)),min Tem — }, 



□ 
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1.9 solitons are limit at infinity 

Definition 1.9.1 A solution g(t) to the normalized Ricci- Hamilton flow is called a 
Ricci- soliton if there exists a one-parameter family of diffeomorphism (f>(t) : M — > M 
such that 

g(t) = <f>(t)*g(0). 

Let X = 4:(j>(t). If X = — V/ for some smooth function, then we call g(t) the 
gradient soliton. 

Clearly we have 

d t g = L x g. 

Then we have we have 

Lxg = (r - R)g. 

This means that X is a conformal vector field. Hence <fi(i) are conformal diffeomor- 
phisms. Note that 

{L x g)ij = V i X j + VjX h 

we get 

(r - R) gij = ViXj + VjXi. 
If g{t) is a gradient soliton, then X = —V/, and we have 

(r - R) 9ij = -2D tj f. 

Taking the trace of both sides we get 

Af = R-r. 

Clearly in the case when M = S 2 , the conformal group of S 2 gives the gradient 
solitons on S 2 . 
Recall 

Mi^DiDjf - l -{R-r) 9ij 
Then we have the following evolution for Mjj 

Lemma 1.9.2 Along the normalized Ricci- Hamilton flow, we have 

(d t - A)M ij = -2RMij + rMij. 

Proof: We do computation in normal coordinates. Recall that 

Af = R-r 

and 

f t -Af = Rf + b 
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where b is a constant. Then we compute 

1 
2 

and 



Using the relation that 



ihJu = Dy t - (r* ),/*. 



Rijki — -^(gudjk — 9ik9ji) 



we obtain 



DijAf — f kkji — f k j k i — Vi(Rjifi) 

= fkjik — Rikjfki + Rufij — Rjifa — ^iRjifi 
= fjkik + Rikjfki + Rufij — Rjifu — ViRjifi 
= fjikk - ^k(R l jk ifi) - R\kjfkl 
+ Rufij — Rjifu — ViRjifi 

= Afji — -(Rifj + fiRj — RkfkQij) 
2R(f ij -^Afg ij ) 



By the definition of My, we have 



Compute, 



Mij = fij - ^{R-r)gij. 



Wv = d t {Dlf- l -{R-r) 9ij ) 

= D%j t - (r*)*/* - \Rt9i 3 + \{r- R)d t9l3 

= D^Af + rf) + \(V t R5* + VjRS* - V k Rg l3 )f k 

-\{AR + R(R - r))g l3 + \{R~ r) 2 9ij 
= DijAf + ^{ViRVjf + VjRVif - V k RVkfgij) 

-^ARgij + rMij. 



Combining these relations together we get 

1 



dtMij = Afij - -AR 9ij + (r - 2R)M, 



2 avj ' v /-^IJ 
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This implies that 

dtMij = A (jij - 1 -{R - r)g^j + (r - 2R)M ij . 

□ 

From this lemma we can easily get 

Proposition 1.9.3 Along the normalized Ricci- Hamilton flow, we have 

(d t - A)\M tJ \ 2 = -2R\M tJ \ 2 - 2\D k M tJ \ 2 . 

Proof: It is easy to see that 

d t \M l3 \ 2 = 2M ij dM ij 

= 2M ij (AM ij + (r - 2R)M ij ) 

= A\MiA 2 - 2R\M i:j \ 2 - 2\D k M ij \ 2 

□ 

This is our result mentioned in section 1.4. 

By the maximum principle we get the important estimate 

Proposition 1.9.4 Along the normalized Ricci-Hamilton flow, there exist constants 
Ci, Ci > such that 

\M i:j \ 2 < de- C2t . 

Recall that our Ricci-Hamilton flow is a family of the metric g(t) satisfying 

dt9ij = (r - R)gij. 

Set 

9 = $9 

where <pt is the diffeomorphism generated by the gradient vector field V/. 
Then we have 

dtg^ = 2M ij (g) := M tj 
From the estimate above we immediately know that 

\Mij\ 2 < C ie ' C2t . 
Then we have the limit at infinity for Mif 

(Moo)ij = 0. 

So we have 

-Roo = r. 

This means that R converges exponentially to the constant r. From this we have that 
the scalar curvature R converges exponentially to the constant r. This eventually 
implies that the metric g(t) converges exponentially to the constant curvature metric 

9oo- 
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Chapter 2 

Bartz-Struwe-Ye estimate 



We first introduce the Moving plane method. Then we discuss the beautiful argu- 
ment of the new Harnack inequality found in [4] for the Ricci-Hamilton flow on the 
sphere S 2 . In this kind of argument, they have to use the uniformazation theorem 
on the sphere S 2 . Note that in the first Chapter, the uniformazation theorem is a 
by-product of the convergence theorem proved in the last section. 

2.1 Moving plane method 

In this section, we introduce an application of the maximum principle. By now it is 
called the moving plane method. This tool was first introduced by A.D.Alexanderov, 
was used by J.Serrin to study a free boundary problem. More than 20 years ago, it 
has been used by Gidas, Ni, and Nirenberg [16] to prove the following result: 

Theorem 2.1.1 Given a locally Lip schitz function f : R R. Assume that 

u e C 2 (B 1 )nC 1 (B~ 1 ) 

satisfies the following partial differential equation: 

—Au = f(u), u > in Bi 

and 

u = 0, on dB 1 . 

Then u = u{r) forr = \x\, i.e., u is a radially symmetric, so that u'(r) < in (0, 1). 

If n = 1, the theorem above implies that u is a even function in [—1,1]. 

Proof: We will only prove the difficulty case when n > 3. 
Choose a constant K > such that 

\f(y)-f(z)\<K\y-z\ 
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for all y,z G [— M, M] where M = maxgj \u\. Since our equation is rotationally 
invariant, we need only to prove that u is reflection symmetric. We will prove that 
u is symmetric with respect to the plane x 1 = 0. 

We now consider the reflection point x x of the point x e B\ with respect to the 
plane x\ = A. Then for x = (xi,X2, ...,x n ), we have 

x x = (2A - x 1 ,x 2 , -,x n ). 

Let 

S A = {x E B 1 ;x 1 < A}. 

and define 

u x (x) = u(x x ). 

Then we have the equation 

-Au x = f(u x ), in Sa- 
lience we have 

-A(u\ -u) = f(u\) - f(u) := C x (x)(u x - u) 
on £ A , where |C A (a;)| < K. Note that 

u\ — u = 0, on <9E A 
Let v = sup{u — u\, 0}. Then we have 

/ |Vw| 2 = - / Av ■ v 
= - / (\U)r 2 



= K v 2 
Jt, x 

The left side is bounded by 



By the Sobolev inequality we can bound it from above by 



< K\Zx\ 2/n / |Vw| 2 
Jt, x 

So for |£ A | sufficient small, we have the upper bound 
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Hence we have v = on such a Ej. So we have 

u\ > u, on T,\. 
Thus we have proved the following 

Assertion 2.1.2 There exists a e\ > such that for all A G (— 1,— 1 + ei), we /iove 

U\ > u, on T,\. 

Define 

A = sup{fi G (—1, 0);u\ > u; in Sa for all A < //} 
from the assertion above we have 

-1 < A < 0. 

Assertion 2.1.3 

A = 0. 

In fact, if A < 0, then we will show that there is a e > such that for all A G [A, A + e) 
we have 

u\ > u, on T,\. 

This implies that A can not be the supremum. Thus we get a contradiction. 
We now choose a compact set D s C £t for some S > such that 

\Xj-D s \<5. 

Then we have 

uj — u> a(8) > 
on Ds- So by continuity, we can find a e± — ei(S) > such that for all 

A G [A,A + ei) 

it holds 

u\— u > a(S)/2 

on Ds and 

|S A - D s \ < 25. 

Since 

-A(u x - u) = Cx(x)(u x - u) 

on Sa — Dg and 

u\ — u > 

on d(E\ — D s ), we can use the argument for Assertion 4.7.2 to prove that 

u\ — u > 

on S A — Ds for 5 > small enough. Thus we have A = 0. □ 
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We close this part by recalling the Hopf boundary point lemma. 

Theorem 2.1.4 Let VL be a bounded domain in R n . Let u G C 2 (f2) n C 1 (r2) be a 
solution of the following partial differential inequality: 

(d - A)u + c(x,t)u < 

where c(x,t) > is a bounded function on Q x [0,T). Assume that there is a point 
x° G dQ and a time t G (0,T) such that 

u(x°, t ) > u{x, t) 

for all (x,t) G Q x (0,T) and u(x°,to) > 0. Suppose that Q satisfies the interior 
sphere condition, that is, there is an open ball Bcfi such that x° G dB . Then we 
have 

jr (x °,t >0 
ov 

where v is the exterior unit normal of dQ at x°. 
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2.2 A new Harnack estimate 

Recall that we study the Ricci-Hamilton flow on the sphere S 2 , 

9t = (r- R)g, (1) 

with g denotes the evolving metric, R the scalar curvature and r the average of R. 
Let c := g S 2 be the standard metric on S 2 of volume 4n. Let 

9 = e u c. 

Then the flow (1) becomes 

de u 

— = A c u-2 + re u , (2) 

on S 2 where A c is the Laplacian on the standard metric c. 
We consider S 2 as the standard unit sphere in R 3 and let 

F : S 2 - R 2 

be the stereographic projection , whose inverse is given by 

F-\x) = (2x, \x\ 2 - 1)/(1 + |x| 2 ), x e R 2 
We introduce the following coordinates 

x = (X!,X 2 ) 

around the north pole po — (0, 0, 1): 

G : R 2 - S 2 - { Po } 

G(x) = (2x, 1 - |x| 2 )/(l + |x| 2 ), \x\ < 1 
For a given smooth function f on S 2 , we set 

(F- 1 )*(/^ 2 )=7te 2 , 
where g^? denotes the Euclidean metric. Then / is given by 

J(x) = f(F'\x))A/(l + \x\ 2 ) 2 (3) 
A simple computation leads to 
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Lemma 2.2.1 Define 



ao = /(po) = (/•<?) (0, 
a, = d(f-G)/dxi(0) 

and 

= d 2 {f ■G)/dx i dx j (0). 
Then we have the following expansions near oo: 

f(x) = A/lx^ ^a + aiXi/\x\ 2 + {^aij -2a 5i j )x i x i /\x\' i + 0(l/\x\ 3 ) 

Wo* = -^(°° + 7? 

+ W'W a ' x ' + 0{1/M7) 

We will call these expansion as (4). 

Let Mo be a smooth function on S 2 , and letw be the unique smooth solution of 
(2) with initial value uq on a maximal time interval [0, T*). We let / = e u , define 
w, f in terms of 

(F- 1 y(fg s2 )=Jg R2 , J = e w 



and let 



Note that 



ao(t) = f(Po,t), 



4 



= ix(F~ 1 (x, £)) + log - (5) 
and that w satisfies the following flow equation 

e w ^ = Aw + re w (6) 
at 

where r is the average scalar curvature of the metric 

e u g s *- 

By Lemma 2.2.1, the expansion above holds for f(.,t) with ao = ao(t),aj = Oj(t) 
and aij(t). we notice that the expansion is uniform for all x G [0, T], where T is any 
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given number in [0, T*). Our purpose is to estimate Du. We introduce the mass 
center y(t) of w(.,t) as y(t) = (yi(t), y2(t)), where 

yi{t) = 



4a (t)" 

Then we have the following 

Proposition 2.2.2 The is a constant C > depending only on ||-u ||c 4 su °h that 

\W)\\<c 

for allte [0,T*). 

Proof: Given T G (0, T*). Doing a rotation of coordinates and the transformation 
x 2 — > —x 2 if necessary , we may assume y 2 (T) = maXi\yi(T)\. By the expansion (4) 
and the arguments for Lemma 4.2 in [16] we derive that for some Ao > 1 depending 
only on ||-u ||c 4 the following holds: 
For each A > A , 

J(x,0)>J(x x ,0) 

whenever x 2 < A, where x x = (x l7 2X — x 2 ) for x = (^i,^)- (Note that x x is the 
reflection of x about the plane x-i — A.) Hence, 

w(x, 0) > w(x x , 0) (7) 

whenever x 2 < A , A > Ao- 

Using the same argument and the fact that the expansion (4) for /(.,£) is uniform 
for all t G [0, T], we can find some Ai > A such that for each A > A x 

w(x,t) > w(x x ,t) whenevert G [0, T] and x 2 < A. 8) 

We are going to show that y2{T) < Ao- To prove this, we consider the function 
u> A (a;, t) := w(x x , t) on the region rr 2 < A, < t < T and define 

I = {A : A > A ,A > maxQ< t < T y 2 {t),w x < A} 

Note that w A solves the equation (6) and coincides with w along the plane x 2 = A. 
By (8), / is nonempty. 

I is also open. Indeed, w x := w can never happen for A > A because of (5). 
Hence for a given A G /, the maximum principle implies 

w x < A for aj 2 < A, (9) 

and the proof of the parabolic version of the Hopf boundary point lemma implies 
that 

49 



along the plane x 2 = A. Consequently, 



f X <f for x 2 <\ (11) 



and 



along the plane x 2 = A, where 



|^ < (12) 

OXo 



f\x,t) :=f(x\t) 



is defined on x 2 < A. For each fixed t G [0, T], we can move the origin to y(t) and 
find the new expansion for / at time t in the following form: 

= ^(ao + ^^ + 0(M- 3 )), (13) 

and 

16 

d x if = -^a ^ + 0(|a:|- 7 ), (14) 



with different coefficients a^. The plane x 2 = A becomes the plane x 2 = A — y 2 (t) in 
the new gauge. Since A G /, we have A — y 2 {t) > 0. Hence we can argue as in [16] 
to show that there is an e(t) > with the property that: 

if s G (A — e, A + e), then f s (-,t) < f(.,t). Using the facts that A > maxo<t<Ty2(t) 
and the expansion above is uniform for all t G [0,T], we can choose e(t) for all 
* G [0,T]. Clearly 

f X < f 

is equivalent to w x < w . By this we obtain that 

(A-e,A + e) c I 

for some e. Thus we have proved the openness of I. 

W now try to prove that I is closed in [Ao, oo). Let A > Ao be in the closure of /. 
By continuity of w, we have w x < w and A > max < t <Ty2{t) ■ If A = max < t <Ty2(t), 
then we have A = y 2 (t 2 ) for some t G [0, T\. Consider now the point yo(t ) as the 
origin and the corresponding stereographic projection F : S 2 — > R 2 . Define z and 
z x by 

F*(e w g R2 ) = e z g s2 , F*(e wX g R2 ) = e zX g s ,. 

Then z, z x are defined on x [0, T] for a semi-sphere 5+. The functions z and z x 
satisfy the equation (00) and we also know that z x < z and z x and z are the same 
on dS+. Using the Lemma above and the expansion (13-14) above at to, we have at 
the north pole N that 

de z de zX 
dv dv 
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where v is the inward unit normal of dS\. this implies that 

dz dz x 
dv dv 

at the point (N,t ). By the Hopf boundary point lemma we have that z x = z on 
dS\. This gives us that w x = w, which is a contradiction with the property that 
w(x, 0) > w(x x , 0). From this we conclude that whenever A G /, it holds that 

A > max < t < T y 2 {t). 

This shows that I is closed. Hence we have that / = (Ao, oo). This means that 

V2(T) < A . 

Then |y(T)| < C. Since T is arbitrary, we proved what we wanted. □ 

It is also clear that the proposition above implies that \Vu(N,t)\ < C for all 
t G [0, T*). By using a rotation, we can transform any point in the sphere into the 
north pole. Thus the proposition gives us the following uniform gradient estimate 

Theorem 2.2.3 Along the Ricci- Hamilton flow (2), we have a uniform constant 
C > such that 

\Vs*u(x,t)\ < C 

for all x G S 2 . 

Integrating along a great circle connecting the maximum point and minimum 
point of e u<yt \ we get that 

inf e u > Csupe". 

M t M t 

where M t = M x {t}. This is a new Harnack inequality for the Ricci-Hamilton flow 
(2). 

With this Harnack inequality we can prove the global existence and convergence 
at infinity as before. 
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Chapter 3 

Hamilton's another proof on S 



In this part, we describe the new proof of R.Hamilton [19] for the convergence of 
the Ricci-Hamilton flow. Usually this is called a proof via an isoperimetric estimate 
and a singularity analysis. 

We study a metric g = {gij} on the two-sphere S 2 evolving under the Ricci- 
Hamilton flow, which is not normalized and on a surface takes the simple form 

— 

where R is the scalar curvature of g. This flow is different from the normalized Rice- 
Hamilton flow just by a rescaling.The reason for this change is that we are consid- 
ering geometric quantities which are invariant under the delation. So these changes 
does not affect any any results claimed true for the normalized Rice-Hamilton flow. 
We point out that in this section, we call this un-normalized flow the Ricci-Hamilton 
flow, which is the only place that is different from the previous sections. 

We now consider any curve A of length L on S 2 dividing the total area A into 
two connected parts A 1 + A 2 = A, and take the smooth curve A of least length L 
on the round sphere of the same total area A = A dividing it into two connected 
parts A 1 + A 2 = A with A 1 = A 1 and A 2 = A 2 the same as before. We form the 
isoperimetric ratio 

C S (A) = L/L, 

and we let 

C s = mfC s (A) 

A 

taking the infimum over all smooth curves A on S 2 . It is well-know that the shortest 
curve A cutting off areas A\ and A 2 on the round sphere is a circle of latitude, and 
its length L is given by 

4tt _ 1 1 
J?~T 1 + T 2 

In fact, in the stereographic coordinates, we can write the metric on the standard 
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sphere as 




Let 7(0) = (r , 6) be a circle of latitude. Then its length is 




We have 





I 2 ^ A 1 A 2 




Using this formula, we find that the isoperimetric ratio of A can be expressed as 



By definition, on the round sphere we have Cs — 1 , and in any other metric 
we can come as close to this as we wish by taking a very short curve like a small 
circle. Hence Cs < 1 in any metric. When Cs < 1 we will show that the value of 
Cs is attained by Cs(A) for some A which is a single simple closed loop of constant 
curvature. The constant C s is also the Sobolev constant, defined as the best constant 
in the following inequality 



where the infimum is taken over all constants c. In fact equality is attained precisely 
when / has a jump discontinuity along the optimal curve A and is constant on either 
side. For a outline of a proof, see the nice survey article by Yau [37]. 

Theorem 3.0.4 (Main Theorem) The isoperimetric ratio Cs increases under the 
Ricci flow on the two-sphere. 

This estimate gives us a new proof of a theorem by Ben Chow [9] . 

Corollary 3.0.5 ([9]). Under the Ricci flow on the two-sphere, any metric ap- 
proaches positive constant curvature. 





inf / (/ - cfda 
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Proof: As we already met in previous sections, by the Klingerberg lemma [7], the 
injectivity radius can be controlled by the maximum of the curvature and the length 
of the shortest closed geodesic circle. The isoperimetric ratio controls the length of 
a short geodesic circle in terms of the maximum curvature (see [8]). Therefore, if 
we take a sequence of points in the space-time where the curvature is as large as it 
has ever been and dilate so the curvature there is one, we can extract a convergent 
subsequence. 

If the curvature times the area is unbounded, the limit will be complete but 
not compact, and we can arrange that it is an eternal solution to the Ricci flow, 
that is a solution defined for all time t G (—00, +00). Since the scalar curvature 
is bounded below, after dilating the limit will have nonnegative curvature. Since 
it is not flat, the strong maximum principle shows it is strictly positive. By the 
Gromoll-Meyer theorem which says that any n-dimensional complete non-compact 
Riemannian manifold with positive sectional curvature must be diffeomorphic to 
R n , we know that the limit is (R 2 ,ds 2 ). It then follows that the limit is the cigar 
solution 

2 dx 2 + dy 2 

ds = T~i — 2~i — 2- 
1 + x z + y z 

In fact, we can write the limit metric as 

ds 2 = g(x, y)(dx 2 + dy 2 ). 

Let 

x = e u cos u 
y = e u sin v 

for u G (—00, +00) and v G (0, 2n}. Then we have the metric on the cylinder, 

ds 2 = g(u, v) (du 2 + dv 2 ) 

where 

g(x,y) = g(u,v)e~ 2u . 

Since the limit metric is a gradient Ricci soliton (see also the Main theorem in [25]), 
the soliton is moved down by translating in the coordinate u. Then we have a 
smooth function / defined on the cylinder such that the gradient of / is just 
for some constant a G R, which generates the soliton metric. Hence we have 

fu = ag, f v = 0. 

This implies that g = g(u) is a function of u and a ^ 0. Thus g(u)e~ 2u is a smooth 
function of x 2 + y 2 = e 2u and g(u + at) defines the soliton. Note that the scalar 
curvature of ds 2 is 

R=-(^] ■ 
9\9 J u 
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Using the equation of the Ricci-Hamilton flow we get 



ag u 



9u 

9 / a 



Solve this and we obtain that 



Ce Cu 



Since g > 0, we must have 
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Ci ± ae Cu 

Ce Cu 
Ci + ae Cu ' 

After rescaling of the constants, we have 

_ e Cu 
9 ~ l + e Cu ' 

Since g(u)e~ 2u is a function of e 2u we can arrange the constant C such that 

9 ~ 1 + e 2 "' 

and this gives us that 

, 9 <ir 2 + dy 2 
as - 



1 + x 2 + y 2 

But the cigar has isoperimetric ratio zero; hence it cannot occur as a limit of surfaces 
whose isoperimetric ratio is bounded away from zero. 

The other possibility is that the curvature times the area is bounded. In this 
case the limit of the dilations will be compact. Again it will have strictly positive 
curvature. Since the normalized entropy 



E = J R\og(RA)da 



is monotone decreasing for the Ricci flow, it must be constant on the limit flow. 
But then Ben Chow's proof (see section 1.6.2) of the entropy estimate shows that 
a certain integral of a positive expression is zero, which implies that the limit is 
a compact homothetically shrinking solution. But we know this can only be the 
sphere. □ 

We now discuss the existence and smoothness of the optimal curve A. 



Theorem 3.0.6 On any surface with Sobolev constant C s < 1 (the isoperimetric 
ratio on the sphere) the minimum is attained on a smooth curve A. 



56 



Proof: If the length L is sufficiently short, the curve lies in a single coordinate 
patch where the metric is bounded above and below as closely as we like by the 
Euclidean metric 5^, say for any e > 

(l-e)5 ij <gij < (l + e)%. 

Let L be the length of the curve measured in Euclidean metric and let A is the 
enclosed area in the coordinate patch measured too in Euclidean metric. Then we 
have 

Vl - eL < L < Vl + eL . 
If Ai is the enclosed area in the coordinate patch, then 

(l-e)Ao < A 1 < (l + e)A . 

Hence, by the isoperimetric inequality in the plane 

Ll > 4vrA , 

we have 

L 2 > —An A, 
~ 1 + e 

and hence 



Assertion 3.0.7 For every rj > there is a A > and a > such that z/C|(A) < 
1 — T) then L > A and Ai > a and A 2 > a. 

In fact, if Cg < 1, we can approximate it as closely as we wish by Cs(A) < l — i] 
for some A. Then this A has length L > A. Since L is not too small, neither area 
can be too small either, so we can find an a > such that A\ > a and A 2 > a. 

Both Ax and A 2 are no more than the total area A, so they are also bounded 
above. Then L is also bounded above; in fact 

i2 (i + £) /4 '= Cs(A)<1 

and 

1 1 2 

A 1 A 2 A 

so 

L < v / 2ttA 

This restricts the geometry of any curve A with Cg(A) < 1. □ 
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The fact that the lengths are bounded does not tell us a lot. If we parametrize 
the curve by arc length s then 



i , dx l dxi 

L = J V^dTdr* 

so in a coordinate patch where gij is comparable to the Euclidean metric, the coor- 
dinate functions x 1 and x 2 have one derivative with respect to s in L 1 , which makes 
them continuous, but this does not help much. Since the dimension is low, standard 
techniques will apply. But here we follow a new approach of R.Hamilton. 

The idea is to take our approximating curve and improve it by running the curve 
shrinking flow on the surface for a short time. The curve shrinking/shortening flow 
is 

where N is the unit normal defined by normalizing dT/ds and k is the curvature 
of the curve T. The existence of the flow can be proved by expressing the nearby 
curve as a graph of the initial curve (see appendix C). Here we move each point on 
the curve in the normal direction with velocity equal to the geodesic curvature k at 
that point. The enclosed areas A 1 and A 2 evolves at a rate 

— / kds, — — = + I kds 



dt J dt 
while the length evolves at a rate 

dt J 

By the Gauss-Bonnet theorem 

J kds + J Kda = 2tt 

on Ai (or A 2 ) since a circle S* 1 on S 2 encloses a disk D 2 . Now K and Ai (or A 2 ) are 
bounded, so f kds is bounded. The Sobolev constant evolves at a rate 

d 2 dL 1 dA 1 1 dA 2 1 dA 

j t og s{ ) = ~ ~ A~ 2 ^f + 

and since L is bounded above and A 1 and A 2 are bounded below, we have 

^logC 5 (A)<0 

provided 

dL 

— « -1. 
dt 
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Hence, we can find a constant E such that Cs(A) decreases when A moves by the 
curve shrinking flow unless 

J k 2 ds < E. 

We now recall a basic property of curve shortening flow. 

Lemma 3.0.8 For the curve shrinking flow on a closed surface, there exists a con- 
stant C < oo such that if we have a solution for < t < 1/M with k 2 < M with 
any M > 1 then 

Proof: Apply the maximum principle to 

«(§) ,+ * 

and the result follows. One can find more detail in Appendix C. □ 

We can flow the solution until k becomes unbounded, and hence as long as 

sup tk 2 (p, t) < 1 
v 

because all the derivatives also stay bounded by the same sort of argument. Suppose 
the equality above is reached first at some time r < 1 at a point £ where 

rk 2 {i,r) = l. 

[If not, L is bounded and k 2 is bounded and we get J k 2 ds < E as before.] Then by 
the first derivative estimate we get 

rk 2 (p,r) > 1/2 

for all p whose distance to £ is less than cy/r for some c > 0. Then 

J k 2 ds > c/^ff 

for some other c > 0. This shows Cs(A) decreases unless r is not too small, and 
then J k 2 ds < E for some E anyway. 

Lemma 3.0.9 For every closed surface and every r\ > there exists a constant E 
such that if Cs(A°) < 1 — r\ then under the curve shrinking flow A evolves into a 
curve A t with 

C 8 (A*) < C S (A°) 

and A* has 

J k 2 ds < E. 
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Proof: By the previous argument, as long as J k 2 ds > E the Sobolev constant 
Cs(A*) continues to decrease. By the previous lemma, we can decrease the isoperi- 
metric constant and bound the L 2 norm of the curvature of the curve along the 
curve shortening flow. This proves the result. □ 

For use of proving our theorem, we also need an algebra lemma below: 
Lemma 3.0.10 For any positive numbers Li, L 2 , Ai, A 2 , A 3 we have 

Proof: If not, we have 
and 

(L, + L 2 ) 2 {^- + - 1 - ) < L 2 2 {^- + - - - ). 
K A 2 Ax + A 3 K A 3 A 2 + Ai 

Rearranging we get 

A 1 (A 2 + A 3 ) < L\ 



A 2 - 






1 


A 2 






1 , 



A 2 (A 1 + A 3 ) ~ (L l + L 2 y 
and 

A^A. + A,)^ L 2 



< 



A 2 (A ± + A 3 ) ~ (Lx + L 2 )2- 
Adding these two inequality we get 

2A 1 A 3 2L ± L 2 



A 2 (A ± + A 3 ) ~ (Li + L 2 ) 2 ' 
Then we have 

2L X L 2 2A l A 3 

' + , , , x < 



(L 1 + L 2 )2 A 2 (A 1 + A 3) 

So we have to have that L\ = or L 2 = 0,and A\ = or A 3 = 0. A contradiction 
with L 1 ,L 2 ,A 1 ,A 2 , A 3 being positive. 
□ 

Now to finish the proof of the theorem, take a sequence of curves A°) with 

C 5 (A°) - C s 

as j — > oo. Deform each by the curve shrinking flow to a curve Aj with 

C 8 {Aj) < C 5 (AJ) 
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and 



k]ds < E 



where kj is the curvature of Aj. Then we still have 

Cs(\j) - C s 
so the Aj are a much better approximating sequence. 



Since each Aj has 



/ 



k 2 ds < E 



we see that the curves are reasonably well behaved in local coordinate charts. Take 
a chart where the metric gij is comparable to the Euclidean metrics 5ij, say < 
9ij < 2(5^, and where the connection forms are bounded, say |r^| < 2. This happens 
for example in geodesic coordinates or in harmonic coordinates. If s is the arc length 
parameter in the metric and T l and N l the unit tangent and normal vectors then 



dJ_ 
ds 



_ rpl 



and 



ds 2 ' * jk ds ds 

Both T and have bounded Euclidean length, so in the Euclidean metric 



d 2 x l 



<C(|*| + 1). 



Since the length L is bounded, we get a bound 

2 





d 2 x l 


/ 


ds 2 



ds<C 



from the bound on k. 

We try to bound the curvature k c of the curve in the Euclidean metric. Let I be 
the Euclidean length parameter. Then 



ds 



dx % dx j 



lJ ds ds 



and dl/ds is bounded above and below since the metrics are comparable. Differen- 
tiating again 

2 dl dH - 25 ^ 

ds ds 2 13 ds 2 ds 
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and it follows that d 2 l/ds 2 is bounded. Then by the rule for inverse functions, we 
have 

d 2 l ds dl . 2 d 2 s 
ds^dl = ~W dP' 
Hence, d 2 s/dl 2 is bounded also, and we get 

2 

dl < C 





d 2 x % 


/ 


dl 2 



in our coordinate patch. This is just a bound 



k 2 ds < C 

for the curvature and arc length of the Euclidean metric. Since in the place 

d6 

k c — ~r 
ds 

where 9 is the angle of the tangent line, we see that 

02-9i = k c ds 

so the angle 9 is of class C 1//2 with respect to arc lengths. Then for a short enough 
segment we can write the curve as the graph of a function y = f(x) with \dy/dx\ < 1. 
Since 

d 2 y/dx 2 
c = [1 + (dy/dx) 2 } 3 / 2 

we now have 

dx 2 ) dX ^° 

and we see y is of class H 2 , and hence C 1+ 5 with respect to x. 

It is now easy to see that since the Aj are locally uniformly bounded in L\ and 
C 1+ 5 we can extract a subsequence which converges in for p < 2 or C 1+a for 
a < 1/2. The limit A will be a genuine immersed curve. A limit of embedded 
curves may not be embedded, but at least it cannot cross itself; at worst it will be 
self-tangent. The limit will still have bounded norm in L\ or C 1+ 2. 

Moreover the limit A attains the minimal ratio Cs(A) = C$- From this and the 
last lemma above it easily follows that it cannot be tangent to itself, for it would 
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be more efficient to take just one part of the curve or the other. Therefore A is 
embedded. Now the usual first variation argument shows A has constant curvature 

2 Ui M 

and hence A is smooth. 

We are now in the right position to give the proof of the Main theorem, which is 
a straightforward calculation. Start with the optimal curve A at time t where C s = 
Cs(A), and construct the one-parameter family of parallel curves A r at distance r 
from A = A on either side. We take r > when the curve moves from A 1 into A 2 , 
and r < when it moves the other way. We then regard L, A±, A 2 and C,s = Cs(A r ) 
as functions of r, and t also by considering the same curves A r at times t near t. 

First we clearly have 

dA 1 dA 2 
— — = L and — - — = — L 
dr dr 

and by a standard formula 



— — =- / kds 

dr 



/' 



where k is the geodesic curvature of the curve A r . (Of course by a standard varia- 
tional argument k is constant on A, but we do not seem to use this, except to check 
A is smooth.) Thus we also get 

— — A\ = [ kds and -r^A 2 = — [ kds. 
dr 2 J dr 2 J 

Now we use the Gauss-Bonnet formula 

Kda+ I kds = 2% 



J^Kda + J 



where the first integral is over the part with area A ± . Differentiating the Gauss- 
Bonnet formula with respect to r gives 

£L d jkds = -j K ds. 



dr 2 dr 



So much for the space derivatives. 

Now we compute the time derivatives under the Ricci flow. This gives 

^ = - / Kds 
dt 

and 

dA 1 

~dT 
dA 2 



dt 



= -2 j l Kda = -4tt + 2 / kds 
= -2 J 2 Kda = -4tt - 2 J kds 
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and of course dA/dt = —8n. So much for the time derivatives. 
Recall that the Sobolev constant of A r is 

/ 1 1 



and taking logarithms for simplicity, we obtain 

log C| = 2 log L - log Ai - log A 2 + log A - log 47r. 
First we compute 

d 

logC 2 

or 

which gives 



d , 2 dL 1 dAi 1 dA 2 

— log C o = 

dr 5 L dr ^ dr A 2 dr 



d , ^ 2 /■ , , , / 1 I 



log C| = - I kds-L . 
dr s s Lj VAi A 2 



At the infimum A when t = t and r = we must have the first variation equal to 
zero. This gives the useful relation 

Zeds = — | — — 
2 

which we can substitute in future calculations. The second space derivative is 

d 2 _ 2 (PL 2 fdL\ 2 1 d 2 ^i 1 fdAiV 1 d 2 A 2 1 /d,V 2 

°£Cg — — - , 2 T2 1 Tx37 ) 4 ,j„2 



dr 2 to 6 L dr 2 L 2 \dr J A 1 dr 2 A 2 \ dr J A 2 dr 2 A 2 \ dr 
which can be worked out to 

d 2 , ^ 2 2d 2 L r2 / 1 1 V r2 / 1 1 

The time derivative is 

(I, 2 dL 1 dAi 1 dA 2 1 dA 

— log C o = 1 

dt b L dt A x dt A 2 dt A dt 

which can be worked out to 

d n2 _2dL 2 /l 1 \ 2 A\ + A 2 2 

dt [ ° gCs ~Ilt' L \M~A~,) AiA 2 (Ai + A 2 ) ' 

Since we have a "heat" equation 

iL r K ds =* L 



dt J dr 2 ' 
so we can get the formula below. 



64 



Formula 3.0.11 The Sobolev constant Cs satisfies 

dt [ ° gCs ~ dr* l ° gCs + A 1 A 2 (A 1 + A 2 ) [1 CsV 

— 2 

Corollary 3.0.12 If C s < 1 it wi/Z increase. 

2 

Proof: It suffices to show that C s + et increases for all e > no matter how small. 
If not, we can find a time t > when it is no larger than at previous times, hence 

— 2 — 2 

at previous times C s + et was no smaller. Hence at previous times C s was larger 
bye(t-t). 

Now pick the optimal A at time t, and construct the family A r as before. Since 
A is a minimum over all r at t — i , we get 

2 

and hence for C s < 1 

|.ogCl>0 

_ — 2 — 

which means at times t < t, C s was not larger than its value at t by more than 
S(i—t) where 5 is as small as we like for t near t. When 5 < e we get a contradiction, 
proving the Main Theorem. □ 
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Chapter 4 

Perelman's W-functional and its 
applications 

Using the trick used in [30], we can give a third proof of convergence of the Ricci- 
Hamilton flow on surfaces. However, we will not do that here. We prefer to give 
some easy but important introduction to Perelman's work [30]. This work is really 
a ground breaking! Before reading this chapter, we suggest the readers look at our 
introduction of works of R.Hamilton in our appendix A of this book. 
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4.1 Is the Ricci-Hamilton flow a gradient flow? 



Some years ago, many people believed that the Ricci-Hamilton flow is not a gradient 
flow. In [30], G.Perelman proved that it is a gradient flow. Actually he introduced 
two functionals. One is called JF functional. The other is the W- functional. Both 
depend on the existence of the Ricci-Hamilton flow. This means that, we first prove 
the existence of the Ricci-Hamilton flow. Then we solve a non-linear back-ward 
scalar heat equation. The functionals depend on the solutions of the non-linear 
back-ward scalar heat equations. Hence, once we know that it is a gradient flow, 
we can expect that there is any period orbit except the fixed points in the space of 
riemannian metrics mod diffeomorphisms. Here the fixed points are Ricci solitons for 
the Ricci-Hamilton flow. In general, classifying Ricci solitons is a difficulty question. 
Let M n be a closed manifold. 

Let g be a Riemannian metric on M and let dV be the volume form of g. Let 
R be the scalar curvature of g. We follow the idea of G.Perelman [30] introducing a 
new functional T as follows 

F= J(R+\Df\ 2 )e- f dV (4.1.1) 
where / is a smooth function on the manifold M. Set 

dm = e~^dfu,. 

We consider the variation for the metric g and for the function /. Write by 8gij = Vij, 
5f — h the variations respectively. Let v = g^Vij, where (g^) be the inverse matrix 
of (gij). Let Rij be the Ricci tensor of g. 
Recall that 

SR = -Av + DiDjVij - R^Vy, 
5\Df\ 2 = v i W i fD j f + 2(Df,Dh), 

5{e- f dV) = { V --h)e- f dV. 

Let dm = e~^dV . We make dm be fixed. Then we have | — h — on M. 
Then we get 

6 J 7 — / [—Av + DiDjV^ - R^Vij + 

J M 

2(Df,Dh) + (R+\Df\ 2 )( V --h)}dm 

Using the facts that 

/ e'fDiDjVijdV = [ DjDie'fvijdV = [ (fjj - f^dV 
Jm Jm Jm 
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and 



we get 



/ e~ f AvdV= [ Ae~ f vdV = [ e~ f [\Df\ 2 - Af]vdV 
Jm Jm Jm 

[ e-f(Df,Dh)dV= [ e-f[\Df\ 2 -Af]hdV 
Jm Jm 

5 J 7 = - Vij(Rij + DiDjf)dm. 
Jm 



Recall that 

/ dm = 1. 

Jm 

then the evolution equations for the gradient flow of T is 



d t9ij = -2(Rij + ViVjf) 

and 

(d t + A)f = -R. 

Let 

X = Vf 

be the gradient vector field on M. Let <p(t) be the one-parameter family of diffeo- 
morphisms generated by X. Set 

9a = {4>*g)ij, f = f ° <t>- 

Then, one can easily obtain that the metric g and / satisfy 

d t9lJ = -2R tJ , (d t + A)f=\Vf\ 2 -R. 
where all geometric quantities are from the metric g In this way we have 

T t = 2 / \Rij + ViVjf\ 2 dm > 0. 
Jm 

This means that along the gradient flow T is non-decreasing in time t. 
For fixed g, we define 

\(g)=m£F(g,f). 
where the infimum is taken over all smooth function with the property 



/■ 

Jm 



e -fdV = I. 



Using the monotonicity of J 7 , we know that X(g) is non-decreasing in time t. 
This property can be used to show that there is no breathers in some special case. 
For detail we ask the reader look at the remarkable paper of G.Perelman [30]. 
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4.2 W- functional and its property 



In this section we introduce the W^-functional of G.Perelman. Then we use it to 
get a lower bound for the injectivity radius of the flow. This conclusion is the little 
loop lemma of R.Hamilton, which is guessed to be by him. This is also one of the 
most difficulty obstruction to Hamilton's program on Ricci-Hamilton flow. With 
this result, one can see that in a compact manifold of dimension three, we can not 
have blow up narrow necks like the shape S 1 x B 2 . 
Introduce the scale parameter r > 0. We define 

W( 9ij J,t) = [ [r(\Df\ 2 + R) + f-n](Anr)- n / 2 e-fdV 

J M 

where / is restricted to satisfy 

/ (A7vr)- n/2 e- f dV = 1. 
Jm 

Note that W is invariant under the simultaneous scaling of r and gij. 
Let a = St. Using the results from the previous section, we get 

5((47TT)- n / 2 e- f dV) = -(v -2h- —)(A7ir)- n / 2 e~ f dV 

2 T 

Then we obtain that 

5W = [ [a(R + |V/| 2 ) + t(v - 2/i) (A/ - |V/| 2 ) 
-TVij(Rij + Dijf) + h 



l(r(R + |V/| 2 ) + / - n)(v -2h- —}{^r)- n ' 2 e-fdV 
2 T 



Take 



Vij = -2(Rij + D^f, ) 

71 

h = -Af-R+-, 
a = -I. 

Taking the trace of the first equation, and then using the second equation, we have 

n 

V = -2(R + Af) = 2(h--) 

At 

Plugging this into the first variation formula for W, we obtain that 

Wt= f [-a(R+\Vf\ 2 -n(Af+\Vf\ 2 ) + 
Jm 

2r\R tJ + D tj f\ 2 -Af-R+ ^-}{^r)- n / 2 e-fdV 

At 
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Since 

we find that 

Wt 



[ Ae~ f dV = [ (|V/| 2 - Af)e~ f dV, 
Jm Jm 

- [ 2T\R ij + D ij f-^-g ij \ 2 (4TTT)- n / 2 e-fdV 



Define 

f*(9ij,T) = miW(gij,f,r) 
where the infimum is taken over all smooth functions satisfying 



/ (Anr)- n/2 e' f dV = 1. 
Jm 



We remark that for fixed g and r > 0, we can use the Sobolev compactness imbed- 
ding theorem, we can find 

Assertion 4.2.1 The minimization problem inf W is solvable. 
Proof: In fact, let 

= e~f . 

Fix the metric g. Then we have 

W(g, f, t) = (47rr)- n / 2 / [4r| V0| 2 - <p 2 log 2 + (tR - n)<j) 2 ]dV. 

Jm 

with phi satisfying 

(47rr)-" /2 / <f) 2 dV = I. 
Jm 

Without loss of generality we may assume that 

Attt = 1. 

So we need only to prove that the following functional 

J((f))= [ [4r|V0| 2 - (f) 2 log (f) 2 + (rR-n)(f) 2 ]dV 
Jm 

has a minimizer in the subset 

A = {0 e H\M); [ <j) 2 dV = 1} 
Jm 

Note that by Sobolev compactness imbedding, A is weakly closed in H 1 and J is 
weakly lower semi-continuous. So to get a minimizer of HJ over A, we only to to 
prove that J is coercive on A. Let u e A. Then using the simple inequality 

logx < nrc 1 /™, V x > 1, 
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the interpolation inequality, and Holder inequality, we have 



/ u 2 \ogu 2 dV < [ u 2 \ogu 2 dV 

JM J{u>l} 



< / u 2+2 / n dV 

J{u>l} 

< [ u 2+2 / n dV 

JM 

< e [ u 2+4/n dV + C(e) f u 2 dV 

JM JM 

<e([ u^dVf? +C(e) 

JM 

<e([ \Vu\ 2 dV + C(e) 

JM 

Choosing e = r we obtain that 

J{4>) = [ {4T\V(j)\ 2 -<p 2 log <P 2 + (tR- n)<p 2 ]dV 

JM 

> [ [3r|V0| 2 + (tR - n)<p 2 ]dV - C 

JM 

> / 3r\W(f)\ 2 dV + mm(rR-n) -C 
Jm m 

This proves the coercivity of J(-). By the direct method in calculus of variation and 
evenness of J(-), we can find positive minimizer of 

J((f>) — > min 

in the subset 

A= {4> e H\M)\ [ (j) 2 dV = 1} 

JM 

One can also minimize J(-) over the subset 

B = {(f) e H\M)\<I> > 0, I (f) 2 dV = l} 

JM 

Here > means that <f)(x) > on M but <fi is non-trivial. As in our earlier work 
[27], we can prove that the minimizer is in the interior of B. □ 

The regularity of the minimizer can be obtained by standard elliptic theory [15] 
One may see [27] and [34] for more about the variation theory on convex subset in 
uniform Banach space. 

Assertion 4.2.2 For an arbitrary metric g on a closed manifold M . the function 
fi(gij, r) is negative for small r > 0. 
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Proof: We follow the argument of G.Perelman [30]. Let f > be a small number 
such that the Ricci-Hamilton flow exists on the time interval [0, f ] for the initial 
metric g^. Let u = (4iTT)~ n / 2 e~~- f be the solution of the conjugate heat equation, 
starting from a 5-function at t — f, r(t) = f — t. Then we have 

W r = -W t < 

and 

W(g(t),f(t),r(t))^0 
as t — > f . Therefore, by the monotonicity of W, we get 

Kg, r) < W(g(0), /(0), r(0)) = W(g(0), /(0), f ) < 0. 

□ 
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4.3 Uniform inject ivity radius bound 

Given a solution Qij(t) of the Ricci-hamilton flow. We begin with 

Definition 4.3.1 We say that gij(t) is locally collapsing at T if there is a sequence 
of times tk — > T and a sequence metric balls B k = B(p k ,rk) at time tk such that 
r 2 /t k is bounded, 



G.Perelman [30] proved the following 

Theorem 4.3.2 If M is closed and T < oo, then the solution Qijit) of the Ricci- 
Hamilton flow is not collapsing at T . 

Proof: Assume that there is a sequence of collapsing balls Bk = B(p k ,t k ) at times 
tk — > T. Then we claim that 



Indeed, we can take 

fk(x) = -\og(f){r^ l dist tk (x,Pk)) + c k , 

where = 4>(p) equals to one on [0, 1/2], decreasing on [1/2,1], and very close to 
zero on [1, oo), and c k is a constant. Note that using 



\Rm\( gij {t k )) < r 2 k , in B k r~ n Vol{B k ) -> 0. 



— oo. 




we get 



e 



f (Anr 2 k )- n / 2 <l> 2 (r^dist tk (x,p k ))dV < (47rr^)-"/ 2 V '6L{B k ) 



J M 



Clearly, since 



r^ n Vol(B k ) -> 0, we have 



— oo. 



Let A k (s) be the area of the sphere S k := dB(p k ,r k s). Let 




be the average of the scalar curvature over Sk- Then we have 



W(9(tk),fk,rl) = 

fjlMmi + (R k (s) + c k - log/ 2 - n)p(s)]A k (s) 
Jo f(s) 2 A k (s)ds 
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As in the last section, we can find a minimizer / of the functional 

fi[4\f'(s)\*- (log f* + n)f*(s)]A(s) 
£ f(s)A( S )ds 

with boundary condition /(0) = 1 and /(l) = 0. Here A(s) is area of the unit sphere 
in the Poincare hyperbolic space H n (— 1). Let m be the corresponding minimum 
value. We now fix this minimizer function / in the all computation in this argument. 
Just as in proving S.Y.Cheng's eigenvalue comparison theorem (the philosophy is 
that the large the domain (or measure), the smaller the first eigenvalue), we can use 
the curvature bound \Rm\(gij(t k )) < r\ in B k to show that 

^W\s)\ 2 - (log f 2 + n)f%s)]A k (s) ^ 

— = < m. 

f(s)A k (s)ds 

Then we have 

W(g(t k ), f k , rl)<m + n(n - 1) + c k . 
Hence, applying the monotonicity of W, we get 

V(g(0),t k + r 2 k ) < p(g(t k ),r 2 k ) -> -oo. 

Since t k + r\ is bounded, this is impossible. □ 

Definition 4.3.3 We say that a metric g^ is k -non- collapsed on the scale p, if 
every metric ball B of radius r < p, which satisfies \Rm\(x) for every x G B , has 
volume at least nr n . 

Clearly this concept is scaling invariant up to a scale. This means that if a metric 
gij is /t-non-collapsed on the scale p, then a 2 gij is /t-non-collapsed on the scale ap. 
With this concept and the theorem above, we can conclude that 

Theorem 4.3.4 Given a metric g^ on a closed manifold M and T < oo, one can 
find a k = K,(gij,T), such that if we have the solution gij(t) of the Ricci- Hamilton 
flow on [0,T] starting at g^ , then gij(t) is K-non-collapsed on the scale T 1 / 2 for all 

te[o,T]. 

Using the convergence theorem of R.Hamilton (see appendix A below) we can 
conclude that 

Corollary 4.3.5 Fix T < +oo. Let g^, t e [0, T) be a solution of the Ricci- 
Hamilton flow on a closed manifold M . Assume that for some sequence t k — > T, 
p k G M and some constant C we have 

Q k := \Rm\(g i:j (p k ,t k )) -> oo, 
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and 

\Rm\(gij(x,t)) < CQ k , 

whenever t < t k . Then up to a scaling factor Q k of metric gij{t k ) at p k there is 
a subsequence of the metrics converges to a complete ancient solution to the Ricci- 
Hamilton flow, which is K-non-collapsed on all scales for some k > 0. 

This result immediately implies that we have a uniform injectivity radius bound 
for the metric sequence involved in the convergent (see also Theorem 4.7in [8]). 
This implies that the little loop lemma (see [21]) of R.Hamilton is true. According 
to R.Hamilton [21], one can accomplish Hamilton program if one can improve the 
Harnack differential inequality. 

Before closing this section, we give some comments on remaining sections of 
[30]: From the section 7 to section 10, G.Perelman generalized the Li-Yau-Hamilton 
Harnack inequality and the Bishop-Gromov volume comparison theorem for Ricci- 
Hamilton flow; from section 11 to section 12, he concentrates to dimension three 
and tries to classify the ancient solution based on Hamilton's pinching estimate. 
These two sections are really hard and really good. His section 13 is about the 
geometrization conjecture of W.Thurston. 



76 



Chapter 5 



Appendix A: Ricci-Hamilton flow 
on Riemannian manifolds 



5.1 preliminary material 

Let (M n ,g) be a Riemniannian manifold of dimension n. In the local coordinates 
(x l ), the metric can be written as 

g = gijdx % dx° 
In short, we just write g = (g^) or g = g^. Let 

(9 ij ) = (.'/;,) ' 
be the inverse of the matrix (gij). define 

d/j, = \J det{g.ij)dx 

be the induced measure (which is the volume form on oriented manifold) on M. 
Let 

r£- = ^{digij + djga-digij) 

be the Christoffel symbols which give us the Levi-Civita connection and covariant 
derivatives of tensor on M. 

The Riemannian curvature tensor is 

p/ — fl.r 1 —fir 1 -\-V l V p —V 1 V p 

Let 

Rijki = dhkRiji, 

We also call this tensor the Riemannian tensor for an obvious reason. 
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We now define the Ricci tensor: 

Rik = Q^Rijkl 

and the scalar curvature function of M: 

R = g lk Rik- 

The Riemannian tensor Rijki is anti-symmetric in pairs i,j and k,l, and sym- 
metric in their interchange. So we can define the curvature operator on 2-forms: 

Rm : A 2 M -> A 2 M 

and 

RlTl^U^ij Q Q RijpqU rs RijpqUpq 

for the 2-form U = U pq dx p A dx q . We have from the symmetry of R^ki that Rm is 
self-adjoint, that is, 

(Rm{U),V) = (U,Rm(V)). 

for the 2-forms U and V. 

Definition 5.1.1 We say that (M,g) has positive curvature operator if Rm is pos- 
itive definite in A 2 M. 

Note that positive curvature operator is stronger than positive sectional curvature. 
The Riemannian tensor R^m satisfies the first Bianchi identity: 

Rijki + Rikij + Riijk = 
and the second Bianchi identity: 

ViRjMh + V ' jRkuh + VkRijih = 

Taking contraction on the second Bianchi identity, we get the following useful 
formula: 

g^ViRjk = l -V k R, 

which is called the second contracted Bianchi identity. 

In dimension three, we have the following important property that 

RijH = gikRji — guRjk — gjkRu + gjiRik — ^R(gikgji — gugjk)- 

This result implies that the full Riemannian curvature tensor is determined by the 
Ricci tensor. Assume that we diagonalize i?^ at a point so that A < fi < v are the 
eigenvalues. Then the only non-zero components of R^ki are those of the form 

-R1212 = ^(A + 

and those derived from it by permutation. Hence, the positive sectional curvature 
condition is equivalent to the one that each eigenvalue of the Ricci tensor is smaller 
than the sum of the other two, which is in turn equivalent to the following condition: 

Rij < -Rgij. 



78 



5.2 Hamilton's program 



In 1982, R.Hamilton denned the following flow for metrics: 

d t gij = -2Rij. 

It is by now called Ricci-Hamilton flow for the metric family (g(t)). It is clear 
that this flow is invariant under the diffeomorphism group of M. Therefore, the 
symmetry of the initial metric is preserved along the flow. The difficulty one has to 
face is how one can prove its short time existence of the flow. However, R.Hamilton 
[18] proved the local existence of this flow on any compact manifold. The proof is 
hard and more involved, and was simplified later by De Turck [13] 

Example: Let M = S n be the standard sphere in R n+1 . Let c^- be its metric. 
Recall that the Ricci tensor for the metric c^- is (n — l)cy. Let 

_ 2 
9ij P Cij ■ 

Since the Ricci tensor is scale invariant, we have Rij(g) = (n — 1)%-. Then we have 
Hence we have 

p 2 (t) = 1 -2(n- l)t. 



note the spheres shrink to a point at time T 



i 



2(n-l) • 



Definition 5.2.1 A solution (g(t)) of the Ricci-Hamilton flow is called a Ricci soli- 
ton if there exists a one-parameter family of diffeomorphisms ($(£)) such that 

git) = &g(0) 

Clearly the example above is a Ricci soliton. Another famous example of Ricci 
soliton is called cigar soliton on R 2 : 

dx 2 + dy 2 

g 



1 + x 2 + y 2 
The generating vector field is 

d d 

x = x jr + y^r- 

ox ay 

Using the Ricci-Hamilton flow, R.Hamilton proved in [18] the following remark- 
able theorem: 

Theorem 5.2.2 (Hamilton's Theorem) Let M be a compact 3-manifold which 
admit a riemannian metric with strictly positive Ricci- curvature. Then M also 
admit a metric of constant positive curvature. 
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We make a remark here about the metric of positive scalar curvature. It was 
asked by S.T.Yau [37] to find a criterion for existence of a Riemannian metric of 
positive scalar curvature on a compact manifold. It is well-know that there is some 
topological obstruction for such a metric. However, the progress is not large if 
we forget the resolution of the Yamabe problem. With this understanding, people 
know that finding a nice Riemannian metric of positive Ricci curvature on a compact 
manifold is a very difficult task. 

R.Hamilton [18] observed that 

Lemma 5.2.3 Along the Ricci- Hamilton flow, we have 

d t T% = -g k \V i R jl + VjRa - V,i^)- 

and 

<9 t log/j = -R. 

where /x = ^det(gij). 

We define B^i = RipjqR kp i q . Using the Bianchi identities we can compute that 
Lemma 5.2.4 Along the Ricci- Hamilton flow we have 

(d t — A)R ijkl = 2{Bij k i — Biji k + B ikj i — B U j k ). 
We can define the square of Rm as an operator 

Rm 2 : A 2 M -> A 2 M 



such that 

Then we can write 



Rvn {U*)ij RijpqRpqrsUrs- 



(Rm ^ijkl RijpqRpqkl- 

In the space of 2-forms on M, we have the Lie bracket defined by 

[U,V] ij = U ip g pq V qj -K P g pq U qj 

Fix a point x G M. This bracket gives us an isomorphism between A 2 M and 
the Lie algebra so(n). So so(n) has an induced metric which comes from A 2 M. 
Let 4> a = {4>p q } be any orthonormal basis in so(n) and let df denote the structure 
constants: 

[0 a , (f) b ] = c abc <j) c . 

We also write c a bc = c abc . For any linear operator 

L : so{n) — > so{n) 
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we define its Lie square as 

= c ace CbdfL cd L e f. 

Computations give that 

(Rm 2 )ijki = 2(Bijki — Bijik) 

and 

(Rm^)ijki = 2(B ik ji — Bujk). 

From these, we obtain that 

(d t - A)R ijM = (Rm 2 ) ijH + (Rm^ijki 

This formula is conceptual more clear than the one above. Using this formula and 
the maximum principle we can conclude that 

Proposition 5.2.5 On the compact manifold M, the positivity of the curvature 
operator is preserved along the Ricci- Hamilton flow. 

Corollary 5.2.6 Under the Ricci- Hamilton flow, we have 

(d t — A)R jk = 2g pr g qs R pjgk R rs — 2g pq R pj R qk . 

We also have 

Corollary 5.2.7 Along the Ricci- Hamilton flow, we have 

(d t - A)R = 2\Rc\ 2 . 

In dimension three, we have 

Lemma 5.2.8 Along the Ricci- Hamilton flow, it holds 

(d t - A)Rij = -Qij 

where 

Qij := 6Sij — 3RRij — (2S — R 2 )gij 

with 

Sij = R 2 j = Ri k g kl Rij 

and 

S = 
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We remark that the tensor Qij = on M 3 if and only if (M 3 ,g) is a three 
dimensional symmetric Riemannian manifold. 

In the rest of this section we assume that M is a compact manifold. Let 



= / Rdfi/ / . 
Jm Jm 



dfi. 

Jm 

be the average of the scalar curvature. 

We now choose a time normalization factor ip — ip(t) such that the metric g = ipg 
has unit volume, i.e., 

dji = 1. 



Then we choose the new time scale t = f ip(t)dt. It is easy to know that 

Rij = Rij, R = tp~ 1 R, f = ip^r 



Since 



we have 



so 



J dfj, — 1, 



M 

log / dji = —r 

M 



dt 



and 



Then we can compute that 



d , 2 

— logv? = -r. 
dt n 



d _ d d , 2__ - 

m 9ij = dt g « + dt logv = n f ~ g -- 2Rij - 

We will call this flow as normalized Ricci-Hamilton flow. It is easy to see that 
these two flows are the same except that a change of scale in space and a change of 
coordinate in time. Along this normalized flow, the volume is preserved. This flow 
is very useful in the blow up analysis for the Ricci-Hamilton flow. Let Rc° be the 
trace free part of the Ricci tensor, that is 

Rc°j = R^ Rgij 

Diagonalize the matrix R^ at a point, and we see that 

\Rc\ 2 = \Rc°\ 2 + -R 2 . 



n 



R.Hamilton [18] proved that 
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Theorem 5.2.9 In dimension three, the positivity of the Ricci curvature is pre- 
served along Ricci- Hamilton flow, and there exists positive constants S and C such 
that the following pinching estimate is true: 

I Rr° 1 2 

JfL < cr-> 

Using this, he proved the famous result mentioned at the beginning of this section. 
In dimension four, he refined the maximum principle, and then proved a similar 
pinching estimate under the assumption that the initial metric has positive curva- 
ture operator. In general dimensions, he proved that the positivity of the curvature 
operator is preserved along the Ricci-Hamilton flow. R.Hamilton conjectured that 
the global existence and convergence of the Ricci-Hamilton flow holds for positive 
curvature operator. Interestingly H.Chen [11] improved Hamilton's result in dimen- 
sion four where only 2 positive curvature operator is assumed. In dimension large 
than four, the conjecture is open. 

In the study of long time existence of the Ricci-Hamilton flow in any dimension, 
we need to do the singularity analysis by using the blow up technique. Clearly the 
Ricci soliton will play a important role in the global problem about convergent of 
the solution metric. It is generally believed that if we have a global solution of the 
Ricci-Hamilton flow, then it converges at infinity to either a Ricci-soliton or a metric 
with constant curvature. 

In analysis of non-linear partial differential equations, blow-up technique is that 
we choose a point and rescale the measure in space-time. Considering it as the 
origin, we then make the following variable change: 

(x,t) -> (\x, \H) 

such that the origin is not blow up point anymore. 

In geometry, for the Ricci-Hamilton flow defined on the time interval A < t < T 2 
where A < and T 2 > 0, we choose a marking on M. By definition, a marking 
is a choice of a point Q G M which we call the origin, and an orthonormal frame 
q at Q at at time t — with respect to the metric g(0). We regard the collection 
(M,g,Q,q) as an evolving complete marked Riemannian manifold. The blow up 
technique can now be defined. 

Definition 5.2.10 We say that a sequence gk,Qk,Qk) of evolving complete 
marked Riemannian manifolds converges to the evolving complete marked Riemannian 
manifold M = (M, g, Q, q) if there exists a sequence of open sets Uk in M containing 
Q and a sequence of diffeomorphisms Fk of the sets Uk in M to open sets Vk in Mk 
mapping Q to Qk and q to qk, such that any compact subset in M eventually lies in 
all Uk and the pull-backs §k of the metrics gk by the mappings Fk converges to g on 
every compact subset of M x (A,T 2 ) uniformly with all their derivatives. 
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Remark: We remark that using cut-off of vector fields, we can complete each F k 
into a global diffeomorphism. 

Example: If F k (x,t) = F\(x,t) = (\x,X 2 t), then 

{F*g) ij (x,t) = \ 2 g ij (\x,\ 2 t). 

Using the delation invariant of Ricci tensor, it is clear that if £ is a solution of 
Ricci-Hamilton flow, so is F£g. Since at the blow up point, we just require that the 
quantity \Rm\ goes to infinity, we have to develop a comparison tool for the behavior 
of the solution at nearby points. R.Hamilton [23] find such a tool for solutions with 
non-negative curvature operator. This tool is called the Li-Yau-Hamilton Harnack 
inequality. 

R.Hamilton ([22]) proved the following result 

Theorem 5.2.11 . Let (M k , g k ,Q k ,q k ) be a sequence of evolving complete marked 
Riemannian manifolds which are solutions to the Ricci flow. Suppose that (1) the 
absolute value of the Riemannian sectional curvature of the M k at time interval 
A < t < T 2 are uniformly bound above by a uniform constant B < +oo, and (2) 
the injectivity radii of the M k at the origin Q k att — are uniformly bounded below 
by a constant 5 > 0. Then there is a subsequence which converges to an evolving 
complete marked Riemannian manifold M = (M, g,Q,q) which is also a solution of 
the Ricci-Hamilton flow, with the sectional curvature and injectivity bound as M k . 

In section 4 in [30], G.Perelman used the VY-functional get the injectivity radius 
bound required in the Theorem above. 

The limiting metric is an ancient solution in the sense that the flow is defined in 
the time interval (— oo,T) for some finite T. So the classification of these ancient 
solutions is a big deal. 

In [24], R.Hamilton studied the non-singular solution of the Ricci-Hamilton flow 
on a compact manifold. By definition, a non-singular solution of the Ricci-Hamilton 
flow is the one which is the solution of the normalised flow existing for all time 
< t < oo, and the curvature remains bounded 

\Rm\ < C < oo 

for all time with some constant C independent of t. R.Hamilton showed that 

Theorem 5.2.12 Any non-singular solution to the normalised Ricci-Hamilton flow 
on a compact three manifold M satisfies one and only one of the following properties: 
C). the solution collapses (that is, the maximum injectivity radius on M at time 
t goes to zero as t — > oo ( we remark that this case is excluded by G.Perelman [30]); 
or 

P). the solution converges to the metric of constant positive sectional curvature; 

or 
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Z). the solution converges to a metric of zero sectional sectional curvature; or 
H). the solution converges to a metric of constant negative sectional sectional 
curvature; or 

H'). we can find a finite collection of complete non-compact hyperbolic three 
manifolds with finite volume Vi, V n , and for all t beyond some time T < oo we 
can find diffeomorphisms 

4>i{t) :Hi^M 

of these manifolds into the manifold M with the solution so that the pull-back of 
the solution metric g(t) by <f>i(t) converges to the hyperbolic metric as t — > oo; and 
moreover if we call the exceptional part of M those points where either the point is 
not in the image of any (pi, or where it is but the pull-back metric is not as close to 
the hyperbolic metric as we like, we can make the volume of the exceptional part as 
small as we like by taking t large enough; and each Hi is topological essential in the 
sense that each (pi injects itiiHi) into 7Ti(M). 

R.Hamilton's argument is clever. R.Hamilton obtained a better scalar curva- 
ture pinching estimate and played the volume comparison of solutions between nor- 
malised and un-normalised Ricci-Hamilton flow in case C),P),Z),and H). In the case 
H'), He used a special parametrization given by harmonic maps. 

In dimension three, G.Perelman [30] (see sections 11 and 12 in [30]) tried to 
improved R.Hamilton's result above. He can give more delicate analysis including 
the volume comparison and other beautiful local and global monotonicity formulae 
of the Ricci-Hamilton flow. It is clear that G.perelman made a great progress in 
the study of the Rici-Hamilton flow. For the developments of Ricci- Hamilton flow 
before G.Perelman's work, one may see the nice survey of H.D.Cao and B.Chow [6]. 
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5.3 The proof of Hamilton's Theorem 

In this section we give a brief proof of Theorem 5. 2.2. 

Let's recall three famous theorems in Riemannian geometry. 

Theorem 5.3.1 (Klingerberg's Theorem) Let M be an even-dimensional, com- 
pact, simply connected Riemannian manifold with strictly positive sectional curvature 
Km- Then the injectivity radius inj on M is bounded from below by n/y/maxM Km- 

Theorem 5.3.2 (Meyers' Theorem) Let(M n ,g) be a complete Riemannian man- 
ifold with its Ricci curvature bounded below, i.e., 

Ric{g) >(n- l)c 2 g > 

Then we have the following diameter bound: 

diam(g) < diam(S n (c~ 1 )). 

Theorem 5.3.3 (Bishop-Gunther-Gromov Comparison Theorem) Let(M n ,g) 
be a complete Riemannian manifold with its Ricci curvature bounded below, i.e., 

Ric(g) >(n- l)c 2 g > 0. 

Then 

V(B x (r))/V c (r) 

is decreasing in r, where V c (r) is the volume function in the space form of constant 
sectional curvature c. 

We now assume that n = 3. Let (e^) be a local moving frame on M 3 . Using the 
volume element form (or Hodge 8 operator) we can define an isomorphism between 
/\ l M and A 2 M such that 

e 1 
e 2 
e 3 

form an orthonormal basis in A 2 M locally. 
Write 

Rm = M pq 9 p A q . 
Using the isomorphism above can choose (ej) such that 

mi 

Rm = (M pq ) = [ m 2 

m 3 

80 























vi" 3 = 


Tf 1 



Then 

1 1 

#1212 = (-Rm(ei,e 2 )ei,e 2 ) = -M 33 = -m 3 , 

etc. By this we have the following expression of the Ricci tensor 

1 / m 2 + m 3 
cm — - I mi+m 3 

y mi + m 2 

The scalar curvature R can be expressed by 

R — mi + m 2 + m 3 . 

By direct computation we find 

m 2 m 3 
mim 3 
mim 2 

Hence the evolution equation for curvature operator c 

(d t - A)Rm = Rm 2 + Rrn^ 

an be written as 

/ mi \ / m? \ / m 2 m 3 
(9 t - A) | m 2 = m\ + mim 3 

\ m 3 ) \ m§ / \ m im2 

According Hamilton's maximum principle we need only to analyze the following 
ODE: 

mi \ / ml 
m 2 I = I ml 




m 3 
We can assume that mi > m 2 > m 3 such that 




and 



mi = sup < MX, X > 
l*l=i 



m 3 = inf <MI,I> 



where M = i?m. Note that m« are globally defined functions. We will consider m« 
as functions of the matrix M. From these definitions we know that mi is a convex 
function of the matrix M and m 3 is a concave function of M. Since 

R — mi + m 2 + m 3 = tr g M 

is a linear function of M.(Here we assume that the metric g is fixed). Hence m 2 + 
m 3 = R — mi is a concave function of M and mi + m 2 = R — m 3 is a convex function 
of M. We have the following results. 
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Theorem 5.3.4 For e G [0, |] ; the pinching condition R > and Rc > eRg is 
preserved along the Ricci-Hamilton flow in dimension three. 

The proof of this theorem is a easy application of Hamilton's maximum principle, 
so we omit it. 

Theorem 5.3.5 For any (3 > and B < +oo ; we can find uniform positive con- 
stants A < +oo and 5 G (0, 1) such that pinching condition j3g < Rc < Bg and 

m l - m 3 < A(m 2 + m^) 1 ^ 6 

is preserved along the Ricci-Hamilton flow in dimension three. 

Proof: By Hamilton's maximum principle, the condition f3g < Rc < Bg is pre- 
served along the Ricci-Hamilton flow in dimension three. 
We now choose 5 G (0, 1) such that 

5mi < m 2 + m 3 

at t — 0. note that this condition defines a convex, closed subset, so we can use 
Hamilton's maximum principle to conclude that 

Smi < m 2 + m 3 

is true for all t > 0. Choose A > such that: 

mi - m z < A{m 2 + m^) 1 ' 5 , (*) 

at t — 0. This is possible since 

mi — m 3 < mi + m 2 < B 

and 

m 2 + m 3 > j3 > 

at t = 0. 

We claim that (*) is preserved along the Ricci-Hamilton flow. In fact, we can 
see that (*) defines a closed convex subset of matrices with 

m i > m 2 > m 3 

and 

m 2 + m 3 > 0. 

So we only need to check that the inequality is preserved by the ODE. 



Compute, 



d m\ — m\ + m 2 m 3 — mim 2 

— log(mi-m 3 ) = 

at mi — m 3 

= m 1 — m 2 + m 3 
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while, 



Hence 



d m\ — m\ + mim 3 + m\m2 

— log(m 2 + m 3 ) = ■ 

at m 2 + m 3 

mi(m 2 + m 3 ) 

> = mi 

m 2 + m 3 



4 log ~r^, — < 5m i - ( m 2 + m 3 ) < 0. 
at (m 2 + m 3 ) L - d 



Therefor (*) is preserved by the ODE. □ 

Obviously we can deduce the Theorem 5.2.9 from the relation 

^ / m 2 + m 3 -mi 

Rc° = - I mi + m 3 — m 2 

y mi + m 2 — m 3 

Using the pinching estimate, we know that the blow time of \Rm\, \Rc\, and R are 
the same. We write the blow up time as T > 0. 

In the following we consider the normalized Ricci-Hamilton flow and we show 
that the normalized metric g(t) converges to a metric of constant sectional curvature. 
Clearly, we need only to show that the normalized metric g(t) converges to a metric 
g with zero trace free Ricci tensor: 

Rc°(g) = 0. 



Lemma 5.3.6 Let g = 4>g with <p > 0. Then we have 



!)■ 


■pfc 

ij 


-pk 
ij' 


2). 


hi _ 

ijk ~ 


" n ijk- 


3). 


Rijkl 


= Rijkl 


4). 


Rij = 


Rij ■ 


5). 


R = (, 


b- l R. 


6). 


dfi = 


(j) n / 2 dfi 



Choose 4> — 4>{t) suc h that 

/ djj, — 1. 
Jm 

Then we have 




and 

n d(j> 
20~di ~ 
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where 

p = I Rdfi. 



M 



Put 



= / <t>(r)dr. 
Jo 



t 

'0 

Then we have the normalized Ricci-Hamilton flow: 

2 

9i9 = -2Rc + -pg. 

n 



Assertion 5.3.7 

fT 

Rmaxdt = +00. 





Proof: Define p by the flow 

d -l -OR „ 

at 

with initial data 

p(0) = i? ma ,(0). 

Then we have 

(d t -A)(R-p) = 2(\Rc\ 2 - R max p) 

< 2(R 2 — RmaxP) 

< 2R max (R — p) 

By the maximum principle we get 

R < P 

on M. Hence we have 

P > Rmax -> +00 

as t — > T. Since for any r < T, we have 

log 44= / 4-!ogP = 2 / 
8 p(0) 7 dt ^ J 

Hence, sending r — > T, we have 



Rmaxdt = +00. 

./o 

□ 



/ 
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Assertion 5.3.8 There exists a uniform constant C > such that 

Rmax — C. 

Proof: Let 

L(t) = diam(g(t)). 

Note that 

V(t) :=vol(g(t)) = l. 
Since Rc > 0, by Bishop-Gunther-Gromov comparison theorem we have 

1 < \.L*. 

On the other hand we have 

Rc = Rc> 2p 2 R mm g = 2(3 2 R mm g. 
So by Meyers' Theorem we have 

f 7T 

L < 



ft V 1 R-min 

By the pinching estimate we have 

Rmax Rmax ^ ^ 

R-min Rmin 

and hence, there is a uniform constant C > 1 such that 

Rmax /-i— 1 



Rmin 



> c- v . 



Therefore, we have 



□ 



7j- 

-Rmaa; ^ C R m i n < C(— y) < C. 



Using this estimate, we can bound V a Rm for any a. Then we can can bound the 
local derivatives of the metric g. Hence we can assume that 

9 -»■ 9, 

and 

R -> c> 

a s t — > T. 
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Note that 

pT pr 

/ r(t)dt= / f(i)t 
Jo Jo 

where 

r = 4>(t)dt. 
Jo 

Since 

and by pinching estimate again, we have 

f f(t)i -»■ +00. 

as r — > T. Hence 

T = +00, 

and by pinching estimate we have 

Rc(g) = 0. 

Actually we can show the exponential convergence. We refer the reader to [17] for 
this. 
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Chapter 6 

Appendix B: the maximum 
principles 

The classical maximum principle is for the scalar heat equation. It says that for any 
smooth solution of the heat flow, point-wise control for at t = is preserved along 
the flow. Assume that M is a compact manifold. A well-known result is 

Theorem 6.0.9 The maximum principle 1. Assume that u is a smooth (C 2 ) solu- 
tion of the semi-linear heat equation 

(d t - A)u = (V, Du) 

where A is the Laplacian with respect to a time- dependent Riemannian metric G(t), 
and V is a smooth vector field on M x [0, T). // we assume \u\ < C at time t = 0, 
then we have 

\u\ <C, Vf > 

A little more advanced maximum principle is the following 

Theorem 6.0.10 The maximum principle 2. Assume that 

u : M x [0, T) — > M 

is a smooth (C 2 ) solution of the semi-linear heat equation 

(d t -A)u=(V,Du)+f(u) 

where D and V are as before, and f is a smooth function on R. If we assume 
C\u < C2 at time t — 0, then we have 

01 < U < 02, Vt > 

where <pi for i — 1, 2 are smooth solutions of the ode 

4> it = f{4>i), 0,(o) = d. 

The maximum principles can be generalized to complete non-compact Riemannian 
manifolds with some restriction (see [1]. 
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6.1 Hamilton's maximum principle 

If f(t) is a Lipschitz function of time variable t. We define 

df , 

~r < c 
dt ~ 

if and only if 

,. f(t + h)-f(t) 
lim sup < c. 

Let F = F(t,y) is a smooth function of the time variable and space variable 
y E R k . Let Y be a compact subset of R k . Let 

f(t) = su P {F(t, y y,yeY}. 

Then f)t) is a Lipschitz function. Then we have the following 

Lemma 6.1.1 (Basic Lemma) 

f t ^4 t F(t,y);yeY(t)}. 

where Y(t) = {y;F(t,y) = f(t)}. 

Proof: Choose a sequence of tj — > t + so that 

lim /feww 

*./->*+ tj — t 

equals the lim sup. Since Y is compact, we can choose yj e V with 

f{t j ) = F{t j ,y j ). 
By passing to a subsequence yj we can assume that 

Vj -»■ J/- 

By continuity we have 

f(t)=F(t,y), 

this implies that y G Hence 

F(t, yj )<F(t,y). 

By this we have 

By the mean value theorem we can find between tj and t such that 

d 

F(tj,yj) - F(t,yj) = ^g(r j ,y j )( y t j - t) 
Sending tj — > t we have Tj — > i and then 



/(t,) - /(t) <| 9 («.»). 



tj-t ~ dr 



This gives us the conclusion. □ 
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Let M be a compact manifold with a Riemannian metric g , and let 



/ = {f } : M —> R k 
be a vector-valued function on M . Let U be an open subset of R k and let 

: U C R k -> R k 

be a smooth vector field on U. We let 5 and depend on time also. Then we 
consider the nonlinear heat equation (PDE): 

|=A/ + W) 

with initial data / = f at t=0, and we suppose it has a solution for some time 
interval < t < T. We let X be a closed convex subset of mappings / : U — > i? fc 
containing the initial data /o, and we ask when the solution remains in A\ To answer 
this we study the ordinary differential equation (ODE): 

on [/, and ask when its solutions remains in X . We define the tangent cone TfX to 
closed convex set X at a point / as the smallest closed convex cone with vertex at / 
which contains X. Then it is the intersection of all the closed half-spaces containing 
X with / on the boundary of the half-space. 

We say that a linear function I on R k is a support function for X at / e <9Af and 
write / G S)# if 1(f) > l(k) for all other k E X. Then <j>(f) £ T f X if and only if 
K<P(f)) < for all / E S f X. 

Lemma 6.1.2 . The solutions of the ODE 

which start in the closed convex set X will remain in X if and only if <p(f) e TfX 
for all f E X. 

Proof: Suppose l(<f>(f)) > for some le SfX. Then 

!ffl = i (| ) = W) )>0 

So /(/) is increasing and / cannot remain in X. To see the converse, first note that 
we may assume X is compact. This is because we can modify the vector field <f>(f) 
by multiplying by a cutoff function which is everywhere nonnegative, equals one on 
a large ball , and equals zero on a larger ball. The paths of solutions are unchanged 
inside the first ball. Then we can intersect X with the second ball to make X convex 
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and compact . If there were a counterexample before the modification there would 
still be one afterward. Let s(f) be the distance from / to X, with s(f) — if / G X. 
Then 

s(f) = supl(f - k), 

where the sup is over all k G dX and all I G SkX. This defines a compact subset Y 
of R k x R k . Hence by the basic lemma 

^ < sup{l{4>{f)) 
where the sup is over all pairs (k, I) with k G dX, I G SkX, and 

s(f)=l(f-k). 

This can happen only when k is the unique closest point in X to / and / is the linear 
function of length 1 with gradient in the direction f — k. 

We now use the fact that is smooth on the compact subset X to get the bound 

Mf) - m\ < w - k\ 

for some constant and for all / and k in X. Then since we have 

l((f)(k)) < 

by the assumption and 

I/- *! = *(/) 

where 

Since s(f) — at t — 0, it must remain for all time. □ 

The key observation of R.Hamilton [17] is that the ODE controls the parabolic 
partial differential equation (PDE). 

Theorem 6.1.3 If the solution of the ODE stays in X , then so does the solution 
of the PDE. 

Proof: As before, we can take X as a compact sunset. Let s(f) be the distance of 
/ G R k from X and let 

s(i) = sup s(f(x, t)) = sup l(f(x, t) - k) 

where the later sup is taken over all x G M, all k G dX and all I G Sk(X). Since 
the sup is taken over a compact set, we can use the basic lemma to get 

j^sit) < sup^l(f(x,t) - k), 
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where the sup is taken over all x, k, I as above with 



l(f(x,t)-k)=8(t). 

Since x is the maximum point of 1(f), we have 

l(Af) = Al(f) < 0. 

Now we have 

j t l(fM — k) — l(Af) + IW)) < TO)- 
Since l(4>(k)) < by assumption, we have some constant C > so that 

iw)) < to) - to) < i0(/) - m\ <c\f-k\ = cs(t). 

Thus we have 

f t s(t) < CS (0- 

Since s(t) = at t = 0, we have s(t) for all t > 0. This shows that f(x,t) remains 
in □ 

We can generalize this result to vector bundles. Let V be a vector bundle over 
the compact manifold M, and suppose V has a fixed Riemannian metric h. Let g 
be the Riemannian metric on M, and let A be a connection on V compatible with 
h. Both g and A may depend on time variable t. We can form the Laplacian of 
a section / of V as the trace of the second covariant derivative with respect to g , 
using the connection A on V and the Levi-Civita connection T on TM for g. Let 
W be an open subset of V and let be a vector field on W tangent to the fibers . 
Then we consider the following nonlinear heat equation (PDE): 

|=A/ + W) 

Let X be a closed subset of U C W. We ask when solutions of the PDE which start 
in X will stay in X . We need to impose the conditions that X is invariant under 
parallel translation by the connection A at each time , and that each fiber of X is 
convex. Then we can judge the behavior of the PDE by comparing to the following 
ordinary differential equation (ODE): 

in each fiber. 

Theorem 6.1.4 If the solutions of the ODE's in each fiber remain in X , then the 
solutions of the PDE remain in X too. 
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Proof: Again modifying the equation we can assume that X is compact . Using 
the metric h in the fiber and writing \ f — k\ for the Euclidean distance from f to k 
in the metric h , we let s(t) be the maximum distance of any f(x, t) from the set X. 
Then 

s(t) = supl(f(x,t) — k), 

where the sup is taken over all x G M, all k G dX in the fiber over x , and all 
support functions / G SkX at k in the fiber at x , The set of all such pairs (k, I) is a 
compact subset of V © V*. Then as before, we have 

^Ksup^lUix^-k) 

where the sup is taken over all x where the distance in the fiber from f(x, t) to x is 
maximal, k is the unique closest point in X to f(x,t) , and I is the linear function 
of length 1 on the fiber of V at x with gradient in the direction from k to f(x,t). 
Again 

^Z(/(s,f)-*) = i(A/) + JM/)) 
and since l((f>(k)) < by assumption 

iW)) < 10(7) - m\ <c\f-k\ = Cs(t), 

where C is some constant bounding the first derivative of on a neighborhood of 
X. If we extend a vector in a bundle from a point x by parallel translation along 
geodesies emanating radially out of x, we get a smooth section of the bundle such 
that all the symmetrized covariant derivatives at x are zero .We extend k G V and 
I G V* in this manner . Since the metric in V is invariant we continue to have \l\ = 1, 
and since X is invariant under parallel translation we continue to have k G dX and 
I a support function for X at k . Therefore 

l(f(x,t)-k)<s(t) 

in the neighborhood . It follows that l(f(x,t) — k) has its maximum at x ,so 

Al(f(x,t)-k) < 

at x . But k and 1 have all their symmetrized covariant derivatives equal to zero at 
x , so /(A/) < at x . From this at x we get that 

^(/). 

Since s(f) — at t — 0, it must remain for all time. This completes the proof . □ 
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In our applications to the Ricci-Hamilton flow, we have a principal G-bundle P 
over M where G is a compact Lie group , E is a vector space with a metric and G 
acts on E preserving the metric , is a G-invariant vector field on E, and Z is a 
closed convex subset of E invariant under . Then solutions of the equation 

|=A/ + W) 

for sections / in P Xq E remain in the set X = P Z. 

Another important observation is that R.Hamilton [20] generalized them to par- 
abolic systems for tensors. As in the matrix case, we introduce the following concept. 

Definition 6.1.5 We say that a symmetric tensor M^- > if MijV l v^ > for all 
vectors vK 

As usual, we let 

A My = g p<1 V p V q M tJ 
be the Laplacian of the tensor on the Riemannian manifold (M,g). 

Definition 6.1.6 Assume that = piMi^gif) is a polynomial in formed by 
contracting products of My with itself using the metric gij . We say that satisfies 
the null-eigenvector condition if whenever is a null-vector of M^-, i.e., 

M ijV l = 0, 

for all j , then we have 

N ij v i v j > 0. 
R.Hamilton [17] established the following 

Theorem 6.1.7 Let (M,g) be a compact Riemannian manifold. Let u l be a vector 
field on M. Suppose that on [0, T], 

(d t - A)Mij = u k V k M %] + Nij, 

where = p(Mij,gij) satisfies the null-eigenvalue condition above. If My > at 
t = 0, then it remains so on [0,T]. 

For a proof of this result, one may consult [17]. Using this maximum principle, 
R.Hamilton [17] proved that in dimension three, non-negativity of the Ricci tensor 
( and positivity of sectional curvature) is preserved along the Ricci-Hamilton flow. 
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Chapter 7 



Appendix C: Curve shortening 
flow on manifolds 



We first derive the standard first and second variations of arc-length functional. 
Then we study the curve shrinking/ shortening problem in a Riemannian manifold. 

7.1 first and second variations of arc- length 

Let (M n ,g) be a Riemannnian manifold. Let 7 : [a, b] — * M be a smooth curve 
parametrized proportional to arc-length, i.e., |7'| = I. 

A variation of the smooth curve 7 : [a, b] — > M is a smooth mapping: 



F : [a, b] x (-e, e) —> M 



such that F(s, 0) = 7(s). 
Let 



X 



F*D U 



Y 



F*D 2 , 



where D 1 — 4- and D 2 



_d_ 

d/i 



. Let 



!,{■)= F{-^). 



We call Y the variation field of 7 
Then we have 
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Hence, at /i — 0, we have 

1 

7 



This formula is called the first variation of arc-length. The critical points of arc- 
length is geodesies in M. 

Assume that 7 is a geodesic. Consider a two-parameter variation of 7: 

F : [a, b] x (-e,e) 2 -> M 
such that F(s,0,0) = 7(5). Let 

X = F.D U Y 1 = F*D 2 , Y 2 = F*D 3 . 



and let 

Then we have 



9 L(^) = f]Y\ {VxY2 ' X)ds 



dfx 2 



and 
d d 



J*W.) = f ^-((V Yl V x Y 2 ,X) + {V x Y 2 ,V Yl X))ds 

J CI I 



I T^ ( yxY 2 ,X){V Yl X,X)ds 

J CI I I 

f J ^-((R(Y 1 ,X)Y 2 ,X) + (V x V Yl Y 2 ,X) + (V x Y 2 ,V x Y 1 ))ds 

Ja |^ I 

/ rij3< v ^,x)(v^i,x)d s 



Thus at (/i!,/i 2 ) = (0,0), we have 



<9/ii9/i 2 v ™ y Z 
+ 



/ (V X Y 2 , V X Y ± ) - (R(Y l , X)X, Y 2 ) 

J a 



This is the second variation formula for the arc-length. 
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7.2 Curve shortening flow in Riemannian mani- 
folds 

We collect some basic facts about curve shortening flow in Riemannian manifolds 
. We believe a good understanding for this flow is also helpful in studying Ricci- 
Hamilton flow on manifolds. One may consult [14] and [12] for more material. 
On an n-dimensional Riemannian manifold (M,g), let 

7 : S 1 x (a, b) -> M 

be an evolving immersed curve. Denote by j t the associated trajectoryi.e. 

7t(0 =-li;t). 

Then the length of % is 

Hit) = [ \-^-^ft\du= [ \^-\du = [ vdu, 
J si du J si du J s i 

where 

du 

is the speed of the curve 7 t . We can define the arc length parameter s by 

d_ ld_ 
ds v du 

which implies 

ds = vdu. 

As usual, we denote by T the unit tangent vector field of 7*, i.e. 

T = . dj = ldj 

ds v du 



Then the time derivative of length is 

d , „ s f dv 

di 



m = Ltt du 



s 1 

d'y 

<V tl ^,T>du 
s i du 

<V u ^,T>du 
dt 



[ < ^,V u T>du 

f < °1 > ds 

7si dt' ds 
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If we require 7 to evolve according to the equation 

97 DT 



then we find that 



where 



dt ds ' 
|i( 7 «) = - / k 2 ds < 0, 



e=-.' DT - 



ds 

is the curvature squared. This leads us to give the following 

Definition 7.2.1 For a curve shortening flow, we mean an evolving immersed curve 
7(-,t) satisfying the evolution equation 

<9 7 _ DT 
dt ds 

Obviously, we can regard 7t(5' 1 ) as a 1-dimensional sub-manifold of M. With 
the induced metric from M, its mean curvature vector field is 

H = (V T T) ± . 

Note that < T, T >= 1. So we have < VyT, T >= 0, which implies VtT _L T^ t . 
Therefore H = VtT. This shows that a curve shortening flow is a mean curvature 
flow of a one dimensional sub-manifold in M. 

Next we will give some fundamental formulae for the curve shortening flow. 

Lemma 7.2.2 The evolution of v is 

— ~ -k 2 v 
dt~ 

Proof: By definition, 

2 (97 (97 
du du 

Differentiating it with respect to t,we get 

dv „ (97 (97 

— = 2 < V< — , — > 
dt du du 

- 2 <V d l d l > 

~ z<Vu dVdu > 



= 2v 2 <V T — ,T> 



DT 

ds 



2 ^ 2 ^ DT DT 

ds ' 9s 

-2A;V 



□ 
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We often need to exchange the order of derivatives. We have 



Lemma 7.2.3 Covariant differentiation with respect to s and t are related by the 
equation 

DT 

VtV s = V s Vt + k 2 V s + R(T, — ), 
where R is the curvature operator on M. 



Proof: We compute 



and 



so 



v s = -v tt , 

V 



d 1 1 

V t V s = — (-)V tt + -v t v u 
at v v 

= k 2 -V u + -V u V t + -R{^^) 
v v v au at 

DT 

= V s V t + k 2 V s + R(T,—). 



Here we have used Lemma 7.2.2. □ 



Lemma 7.2.4 The covariant differentiation of the unit tangent vector to the curve 
with respect to time t is 

D 2 T 

V t T = k 2 T + 



ds 2 



Proof: The proof is a straightforward calculation. 

V«T = V,(^) 
v au 

= k ^ + - V «^7 
v au v at 

, 2^ D 2 T 

□ 
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7.3 Bernstein type estimates 



In this section, we will assume M is a locally symmetric space, i.e. VR = 0. For a 
locally symmetric space,we have 

VR(X,Y,Z,W) = (V X R)(Y,Z,W) 

= V X (R(Y, Z, W)) - R(V X Y, Z, W) - R(Y, V X Z, W) 
-R(Y,Z,V x W)=0 



V X (R(Y, Z, W)) = R(V X Y, Z, W) + R(Y, V X Z, W) + R(Y, Z, V X W), (1) 

for all X,Y,Z,W G TM. We will also assume that M satisfies Condition(X), i.e. 
there exists a positive constant A,s.t. 

R(X,Y,Z,W) < A, 

for all unit vectors X, Y, Z, W. 

With these assumptions, we can give some precise estimates which will bound the 
evolution of 

,D n T, n 



1 ds n 1 

This type of estimate also appears in [32] for the Ricci-Hamilton flow. 
First, let us compute the time derivative of 

D n T 2 
ds n 

as follows: 

-(\ — \ 2 ) 
dt [l ds n 1 ' 

D D n T D n T 

— 2 < > 

dt ds n ' ds n 

DDD n ~ l T ~D n T . DT D n ~ l T D n T 

— o < I- k I- R(T 1 > 

ds dt ds™- 1 ds n y ' ds' ds n ~ l ' ds n 

DDD n - l T D n T 2 D n T 2 DT D n ~ l T D n T 

< d~sdt ds™- 1 , ~ds~^ > + + ^ ' ~ds~' ds™- 1 ' ~£hF> 

D D D D n - 2 T 2 -D n-1 T , DT D n - 2 T D n T 

2 < ds^dsdt ds n - 2 + ds"- 1 + ( ' ~ds~' ds n - 2 >'~d!F > 

, 9l D n T l9 . DT D n ~ x T D n T . 

+ 2fc — — + 2R(T, — — , — r ,^—) 

1 ds n ds ds 11 - 1 ds n 

D 2 D D n - 2 T D n T D,~D n - l TD n T , ~.D n T l9 

— 2 < > +2 < ffc 1 > -\-2k \ 

ds 2 dt ds n ~ 2 ' ds n ds K ds""- 1 h ds n 1 ds n 1 
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D DT D n - 2 T D n T DT D n ~ x T D n T 

_ 2 D n D D n T D i 2 D n -*T D n T 

"' _2< 9s™9t ' ds n > +2 2-^ < di^ ds^ 9s™ > 

i=0 

D i DT D"- 1 "^ D n T 

+ 2 X^ < d~i( R ( T > ds n-l-i )' "^T > • 

Using Lemma 7.2.3, we get 



D n D D n T 

2 < T > 

9s™ dt ' 9s n 

D n+2 T D n T D n T /l9rTlx D n T 

Q s n+2 Q s n Qgn V " 

9 2 D n T 2 D n+l T 2 D n T 2 D n T 



So 



-n— i 2 ^ 

9t U 9s™ 1 J 

a 2 DT 2 L> n+1 T 2 D l 2 D n ^T D n T 

d?^~ds™^ ' 1 9s™ +1 1 + ^ < 9?^ 9s™"* )' > 

i=0 

^ D i DTD 71 ' 1 -^ D n T 

i=0 

It is easy to see that 

£>* ,, 2 Dn ~ iT \_ST C D j+1 T D k+1 T D n -i~ k T 

and 

D* DT D n ~ 1_i T _ ^ . D^T D k+l T D"-i-j-*t 

9^ ( ( ' 97 } as"" 1 -* } ~ ^ CijfcjR( 9^' 9^ r) ds"-i-i~ k ' (3) 

where the coefficients Cj jfc are constants. To obtain (3), we have repeatedly used 
(!)• 

Noting that M satisfies Condition (A), and then putting above equations to- 
gether, we obtain 

d (l D n T ]2) d 2 (l D n T, 2) , D n+1 T 2 D n T 2 DT 2 D n T 2 
dt [l ds n N _ 9s 2 U 9s™ N 1 9s™ +1 1 11 9s" 1 21 9s 11 9s™ 1 
D 2 T D n T 2 ^ D' l T D^T D k T D n T 

0<i,j,fc<n 
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where Cj are positive constants depending on n and A. In the last term, the range 
of indices satisfies in addition either 



i + j + k = n + 2 

or 

i + j + k = n. 
Theorem 7.3.1 Fix t G [0,+oo). Let 

M to =: maxk 2 (-,t ). 

Assume 

M t0 < +00. 

Then there exist constants c\ < +00 independent oft such that for 

fe(f o,io+^log(l + 4Mi ^ A + 1 )], 

,D'T l2< c,M t0 



we have 



ds l 1 - (t-to) 1 ' 1 ' 



Proof: Without loss of generality, we may assume that t — 0, and then translate 
the estimates. 

(1) For I — 1, we have 

d ..DT l9 . D D , 9 L>T , DT. L>T 

a' 1 ^" = 2< aiM T + k l»I + R{T ^ T ^ > 

D 3 T DT 

= 2<—-,— >+4k 4 + 2k 2 R(T,N,T,N) 

OS OS 

8 2 DT D 2 T 

It follows from the maximum principle that M t satisfies 

M t , M 

log -= log -= — < 2At. 

6 |M t + 1 6 I M + 1 " 

If 

^ 2 i 1 ° g(1 + 4M f0 +A + l ) - 

then 

M t < 2M . 



So we may choose c x = 2. 
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So 



(2) For I — 2, we have 

<9 D 2 T 2 2< DDDT k2 D 2 T DT DT D 2 T 

dt ds 2 ds dt ds ds 2 ds ds ' <9s 2 

2 < D 2 DT D 2 T > +2< D ( k 2 DT ) ° 2 > 
ds 2 dt ' ds 2 ds ds ' ds 2 

D . , DT, . T> 2 i 9l D 2 T l9 

T> 2 T 

' (9s ' ds ds 2 

d 2 (] D 2 T 2 \ D3T \2 A c,\ DT \2\ D2T \2 

~ ds 2 ds 2 ds 3 ds ds 2 

D 2 T DT 

o d 2 t , . dt 



:(^l 2 + 3| 



Since 



<9t ds 2 ds 
d 2 /t D 2 T l9 ,DT |9 . ,D 3 T |9 r/ ,DT l9 AX n ,D 2 T l9l 

£ ^('l^l 2 + 31 gjl 2 ) - *l -^l 2 + Wi8| ^l 2 + 2A) - sil^H 

+ (2A* + 12)|f r + 6A|f |». 

i£ 2X log(1+ 4M„ + A A + l )£ 



4M + A + 1 ' ~ 2(4M + A + 1) ' 
we have 

d /t D 2 T l9 ,DT, 9N <9 2 D 2 T i9 ,DT l9N 9 

7w I + 31 9F 1 > 5 ai^ 1 ^ 1 + 31 aF 1 > + 52M » + 12AM °' 

Thus it follows that 

D 2 T DT 

and we may conclude on this time interval that 

,D 2 T |2 < 16M 



1 ds 2 1 - * 
So we may choose £2 = 16. 

The induction hypothesis and repeated usage of the Peter- Paul inequality, i.e., 

ab < ea 2 + —b 2 , 
4e 

allow us to find constants and A, B on our time interval such that 



dt K ^ a% 

i=l 

Thus we obtain c m as before. □ 



i=i 
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Note that these estimates prove the long time existence result. That is, as long 
as the curvature remains bounded on time interval [0, a), one can define a smooth 
limit for the tangent vector T at time a. Thus, by integrating the tangent vector T, 
one can obtain a smooth limit curve. 
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Chapter 8 

Appendix D: Selected topics in 
Nirenberg's problem 

We have two aims in this appendix. One is a density result of the prescribed gaussian 
curvature functions for Nirenberg's problem. This result says that any smooth 
function, which is positive somewhere, on S 2 can be approximated in C 1,cr (for any 
a < 1) norm by a sequence of smooth positive somewhere Morse functions, which 
are the prescribed gaussian curvature functions in Nirenberg's problem. A similar 
result in the scalar curvature problem on the sphere of higher dimension was already 
obtained by Yanyan Li [46], and R.Schoen and D.Zhang [48]. The L p dense result 
in Nirenberg's problem was obtained by J.P.Bourguignon and J.P.Ezin [42]. 

The second result concerns the structure of the solution set of the scalar curvature 
problems (and also Nirenberg's problem) as the prescribed scalar curvature varied in 
a suitable way. The result shows that there are many more free directions at a scalar 
curvature function so that we can find more solutions for the near-by prescribed 
scalar curvature functions. We introduce two methods in studying this problem 
separately. Although these methods are known to many people, it is not apparent 
that one can apply them to our problem here. In fact, to formulate a good condition 
is a difficult and important point. We only formulated a very simple condition to 
scalar curvature problems and to Nirenberg's problem. See our Main Theorem and 
Theorem 3 for the statements. We think that we are the pioneers in such a study 
from the view-point of Catastrophe theory or Bifurcation Theory. Our results in 
this section are new, but we have obtained them many years ago. 

8.1 A C l cT dense result in Nirenberg's problem 

We begin with some basic material. Let (D, ds 2 )be a 2-dimensional manifold. In 
local coordinates (u l ), write 

ds 2 = g = gijdvfdv? . 



Ill 



Let 

\g\ = det(g i:j ). 

For the vector field X = XWj, we define the length of it as 

\X\ = v'.'/oA-'.Y.'. 

We define the gradient operator y of the metric g as follows. For / e C 2 (D), 

let 



Hence, 

lv/l 2 = £ 

Define 



du 1 dui ' 
d 



We define the divergence operator div of g as follows. For the vector field X = 
X j dj, we let 

divX = -^dj(^\X j ). 
V\9\ 

For / G C 2 (D), define 

A/ = diu( V /). 

Hence, 

1 /— •• <9f 



We call A g the Laplacian operator of (D,g). 

In the local moving frame we define the Laplacian operator of (D, g) in the 
following way: 

g = (er + (0 2 ) 2 - 

Recall the structure equations: 

d0 i = i j A6 i ,dBl = -K6 1 AO 2 . 

Define 

du = Uid 1 . 

Taking the exterior differentiation we have 

= dui AP + mP A 6). 



Define Uij by 



UjiO 1 = duj — Ui9 % y 
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By this we have; 

UjS 1 A 9 j = 0. 

Hence Uij = Uji. 
Define 

AgU = U U . 

Obviously one can define A g in higher dimensions. 
In the isothemal corrdinates 

g = exp(2f)(dx 2 + dy 2 ), 

Let 

9 1 = e f dx,9 2 = e f dy. 

It is easy to verify that these two definitions are the same. 

In higher dimensions, we can also verify the equivalent of these two definitions. 

For u G C 2 (D), let g = exp(2u)g. Let/cbe the Gauss curvature of the metric g, 
Let K be the Gauss curvature of the metric g. Then, we have 

Assertion 8.1.1 

K = exp(-2u)(- A g u + k). 
Proof: We use the local moving frame to prove this result. Let 

g = ^) 2 + {6 2 ) 2 . 

Recall the structure equation: 

d& = Q)NB\de\ = —K0 1 A 2 . 

Let 

^ = e u 9\ 

By 

dd* = 6) A & 

we obtain that 

= 9) + Uiffl - Uj0 : . 

So 

d{§\) = d{o\) + d{u 1 e 2 -u 2 e l ). 

By the equation 

d§\ = -K9 1 A 9 2 

we know that 

-Kexp(2u)6 1 A 6 2 = -kO 1 A 6 2 + d Ul A 6 2 - du 2 A 1 + u x de 2 - u 2 d6 l . 
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Note that 



dui = Uijd 3 + UjOj, 



Hence we have 



Kexp(2u)9 l A 9 2 = -kO 1 A 2 + A y u6 l A 9 2 . 



□ 



We now recall the Nirenberg problem in S 2 . Let {S 2 ,c) be the standard unit 
2-sphere in Euclidean space R 3 . Then Nirenberg's problem is to find functions K(x) 
on S* 2 such that they are the Gaussian curvature functions of metrics g, which are 
point-wise conformally equivalent to c on 5 2 . This amounts to finding a smooth 
function u : S 2 —> R satisfying the following equation on S 2 : 



Here A is the Laplacian with respect to c. Assume that u is a solution of (N). 
Then the metric g = e 2u c is conformal to c and has K as its Gaussian curvature. 
Integrating both sides of (N) on S 2 , we find 



So a necessary condition for (N) to be solvable is that K be positive somewhere and 
we write such function class by C^S 2 ). 
We have the following density result. 

Theorem 8.1.2 On S 2 any smooth positive somewhere function can be approxi- 
mated in C 1,a -norm by a sequence of Gaussian curvature functions. 

To begin a proof, let's assume the well-known existence result of A.Chang and 
P.Yang [39], Han [51], and Chang and Liu [52]. The result says that 

Theorem 8.1.3 Assume K is a Morse function on S 2 satisfying AK(x) ^ when- 
ever VK(x) = and K(x) > 0. Let p =number of local maxima and q =number of 
saddle points with AK < 0. If p ^ q + 1, then (N) admits at least one solution. 

So we need only to prove that any smooth positive somewhere Morse function can 
be arbitrarily approximated in C 1 ' CT -norm by a smooth positive somewhere Morse 
function K on S 2 with nonzero Laplacian at each critical points, and with |l+p — q\ 
nonzero. Here p is the number of the critical points of K at which AK < and at 
which the Morse index of K is one, and q is the number of local maximum points 



The construction of such a function can be carried out as follows. Choose any 
minimum point x of a given Morse function H . Then by the standard Morse lemma, 
we know that in some neighborhood of x there is a local coordinate (a; 1 ,^ 2 ) such 



Ait = 1 - Ke 



2u 



(AO- 




of K. 
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that our function H can be written as H(x) = H(0) + (x 1 ) 2 + (x 2 ) 2 . So we can create 
a near-by Morse function H t with only two critical points, one is a local minimum, 
and the other is a local maximum, in this neighborhood, and H t is H outside. This 
H t has one more maximum point than H . Such a construction was also used in [54] in 
the scalar curvature problem in the sphere of higher dimensions (so one may see [54] 
for more detail), and it can be made again at any minimum point of H t . By repairing 
this process we may find an approximated Morse function H t with q as large as we 
want, but with fixed p. Here we measure the approximation in C 1,a norm on S* 2 . We 
remark that the Laplacian of H t at any critical point can be slightly modified to be 
nonzero by enlarging the scale in just one direction, for example, the x 1 -direction. 
Now Chang- Yang-Han-Chang-Liu's existence result above (see also [39] [51] [52]) tells 
us that H t is the prescribed gaussian curvature in Nirenberg's problem. Hence we 
get the C l,a density result for Gaussian curvature functions on S 2 . Such a result 
could not be expected to improve to C 2 norm because Morse functions are stable in 
C 2 function space. 

We point out that our argument for the density result depends on the existence 
result of this geometrical problem. Once we have an existence result for some kind 
of Morse functions, we can use the above construction to get the density result; 
for example, we could prove the density result for the scalar curvature problem in 
dimensions three and four by using the above construction. 



8.2 Bifurcations in the scalar curvature problem 

Suppose (M n ,g ) is a compact Riemannian manifold with positive scalar curvature 
function and n > 3. Given a smooth function K on it. The usual scalar curvature 
problem on M n is to find a conformal metric g to go such that g has the scalar 
curvature function K. We will let iV+1 = (n+2) / (n— 2)+l be the Sobolev exponent. 
Let L = — c n A +R g be the conformal Laplacian so that the first eigenvalue of it is 
positive, where c n = 4(n — l)/(n — 2). Suppose K is the scalar curvature function 
of the metric u^ n ~ 2 go, where Uq is some positive smooth function on M. 

When studying the scalar curvature problem on S n , n > 3, we want to under- 
stand the structure of the solution set. Recall that the scalar curvature problem on 
S n is to find a pair (u, K) e C\(S n ) x C 2 (S n ) such that the function K is the scalar 

4 

curvature function of the metric u n - 2 c, where c is the standard metric on S n . If the 
problem is non- degenerate, then we can try to do the Morse theory (see the famous 
work of A.Bahri and J.M.Coron [40]) or use the standard implicit function theorem. 
We will do some work in this part if it is not in the above case. So we are led to 
the situation where there occurs a "bifurcation" . Following the general direction of 
the lecture note of L.Nirenberg [38], that is, using the Lyapunov-Schmidt process 
to reduce the problem to a finite dimensional one, we try to use the "standard bi- 
furcation theory" and "catastrophe theory" of Thorn-Mather for our problem. The 
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difficulty lies in finding good conditions on the prescribed function such that we can 
use the singularity theory. Our study is motivated by our reading of the works of 
A.Bahri [54], A.Bahri and J.M.Coron [40], A.S.Y. Chang and P.Yang [39] and trying 
to get some solution of small Morse index. 

Suppose (M n ,go) is a compact Riemannian manifold with positive scalar cur- 
vature function and n > 3. We study the scalar curvature problem on M. Our 
main concern in this part is to find some condition on Kq such that there is a class 
of functions near Kq being scalar curvature functions of conformal metrics of g . 
In short, we call such a function a scalar curvature function. It is clear, if K is a 
composition of a conformal diffeomorphism with K Q , then it is a scalar curvature 
function. So, in the sense of trying to use the "catastrophe theory", our problem 
may be considered as finding a suitable condition on K such that K is a composition 
of a conformal diffeomorphism with our Kq. Note in catastrophe/singularity theory, 
they use coordinates substitution to get a normal form. Hence we can not apply the 
catastrophe/singularity theory directly. However, we may first consider this prob- 
lem from the nonlinear functional analysis point of view. Then we will reduce the 
problem to a finite dimension one. 

For simplicity, using a change of conformal metric, we may assume at this mo- 
ment that Mo — 1- 

Let N(u, K) be the operator — Lu + K(p)u N from a second order Holder space 
(C 2 > a (M)) times another Holder space (C°> a (M)) to a Holder space (C°' a (M)). If 
the linearized operator T = N' u (uo, Kq) is a surjective between the tangent spaces ( 
so it is an isomorphism by Fredholm alternativity ), then by the implicit function 
theorem, we know every function near Kq is a scalar curvature function. So we 
assume the kernel of T is a finite dimensional vector space E = sp(vi, ...,v<i), where 
each Vi has unit L 2 -norm. In fact, this is the only case which is of geometrical 
interest. If we introduce the operator F(u,K) = N(u,K) + P(u — uq), where 
P{u — uq) is the L 2 projection of u — Uq to E, then F u (uq, Kq) is an isomorphism and 
we can safely use the implicit function theorem. So we find neighborhoods G and 
G2 of 0, and a neighborhood G\ of Kq in the Holder spaces, and a smooth function 
u = u(f, K) from G x G\ to G 2 such that u(f, K) satisfies: 



Furthermore, our map u(.,K) is a diffeomorphism from Go to G2 for every K G 
G\. Hence, obtaining a solution for the scalar curvature equation is equivalent to 
obtaining a solution of the equation P(u(f,K) — uq) — f on the intersection of E 
with Gq if K is in G\. Now, it seems that for every K e G±, we have a d— family of 
solutions u(f, K), where / is the parameter with range as the intersection of G 2 with 
E, of our scalar curvature problem. So this is the local structure of the solution set 
of the scalar curvature problem if the prescribed scalar curvature function is near Kq 
in the above sense. In fact, this is not true because of the Kazdan- Warner condition. 
For example, if we take 



F(u(f,K),K) = f. 



(Dl) 



M = S n = {x 



(xi, ...,x n+1 e K 



>n+l. 



XI = 



1} 
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be the standard sphere and take K — 1 + ex±, then the scalar curvature problem 
has no solution. 

In the following paragraph, we want to learn more about solving our system (Dl) 
from the view-point of singularity theory. For this purpose, we write 



and u(x, K) = u(^2 XiVi, K). Hence we need only to solve the nonlinear system 



for every i = 1, d. This kind of reduction is well-known, see section two of L.Simon 
[33], for example. We remark that the key to solving this problem lies in under- 
standing the map u(x,K). But it is very difficult to achieve such a understanding. 
Let's go back to solving the nonlinear system. 

So to understand this system, we consider the Taylor expansion of the function 
u(f,K) near (0,Kq). First note that (uo,Ko) is in the zero set of F(u,K). By the 
implicit function theorem we have u(0, K ) = u and u x .(0, K ) = 1 for every i. 

Now, we take the derivatives in the equation (1) and find 





< u(x, K) — Uq, Vi >= Xi. 



Lu x + NKipju"- 1 ^ + P(u x ) = v, 



and 



Lu K + 5K(p)u N + NK{p)u N - 1 u K + P{u K ) 



0. 



From this we find that at (u , K ), 



< 5K(p),Vi > +N < Ku K , Vi > + < u K , Vi > 



0. 



Let's compute the second derivatives and we find that 



Lu xx + NKip)^- 1 ^, + N(N - l)K(p)u N - 2 u 2 x + P(u xx ) = 0, 
-Lu xK + NK(p)u N ~ 1 u xK + NSKu"' 1 ^ 



+N(N - l)K{p)u N - 2 u x u K + P(u xK ) = 0, 



and 



Lu KK + NKip^-'uKK + S 2 K(p)u 



N 



+N5K(p)u N - 1 u K + 



N(N - l)K{p)u N - 2 u 2 K + P{u K k) = 0. 



So at (-u ,-Ko), we have: 



iV(iV - 1) < Ku 2 x , v > + < u xx , v >= 0, 



(D2) 
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iV < SKu x , v > +N(N - 1) < Ku x u K ,v > + < u xK , v >= 0, 

and 

< 5 2 K, v > +N < 5Ku K , v > 
+ N(N - 1) < Ku 2 K ,v > + < u KK ,v >= 0. 

Again taking the third derivatives, we find more relations at (uo,K ). But these 
relations are not useful to us. 

Now going back to our finite dimensional problem and developing the Taylor 
expansion at (uo, K ) we find that 

= 2 < u K , vt > + < u xx , Vi > x 2 + < u xK , Vi > x 
+ < U K K,Vi > +h.o.t.. 

We can consider the right hand as a perturbational system of the first two terms 
because we can let 5K = t 2 h and x = ty with \t\ suitable small. Then we get 

= 2 < u K h, Vi > + < u xx ,Vi > y 2 + t < u xK h, v { > y 

+t 2 < u KK h.h,Vi > +h.o.t.. (-D3) 

So it is a small perturbation of a "quadratic" equation, and by using (D2) we can 
conclude: 

Theorem 8.2.1 Assume d—1. 

(i) Suppose there is a smooth function h such that A :=< Ku 2 x) v >< ux-h,v > 
is nonzero. Then by (D2) and (D3), either A < 0, we have two solutions of the 
scalar curvature problem for such K , or A < 0, we have none. 

(ii) Suppose there is a smooth function h such that < ux-h,v >= but B :=< 
Ku 2 x) v >< u x xh,v > is nonzero, Then by (D2) and (D3) we are in the standard 
Whitney fold case. This means, either B > 0, and we have two solutions of the 
scalar curvature problem for such K, or B < 0, we have none. 

If d — 1 and B = 0, we need to consider the third derivatives. See M.S.Berger 
[54] or Martin Golubitsky and David G.Schaeffer [47] for discussion. We further 
theorize that Theorem 2 (above) is only useful in compact Riemannian manifolds 
which are not conformally diffeomorphic to the standard sphere. 

In the remaining paragraph, we discuss the case when M = S n equipped with 
the standard metric and K = n(n — 1) + H with H is small in the Holder space. 
Here if we take u — 1 again, then we have our d — (n + 2)(n + l)/2 because of the 
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conformal group on S n . So to obtain some interesting result we need some facts from 
the conformal geometry on S n . Now we discuss the conformal Killing vector fields 
on S n first (see [42]. It is well known that any vector w G R n+1 gives a conformal 
Killing field W on S n by 

W(v) = w — (w.v)v 

and this W is the projection of w onto the tangent space of S n at v. Conformal 
Killing vector fields on S n form a Lie group, which is isomorphic to the Lie algebra 
o(n + 1, 1) of the group 0(n+ 1, 1) of all linear isometries of Minkowski (n + 2)-space. 
By definition {vi, v^} is the basis of o(n+l, 1). Given a conformal diffeomorphism 
F : S n -> S n . Then we have F*(u N ' l c) = Tp^^c with T F (u) = \F'\^u o F. 
Note that F*(c) = \F'\ 2 c. This T F is an isometry of the Hilbert space H\S n ). Let 
E be the positive part of the sphere of H 1 (S n ) containing u® — 1. Note the orbit of 

n — 2 

uq = 1 in £ under the conformal group action is {uf '■= \F'\^~}. More explicitly, 
uf are of the forms 

U F (w) = d n (- /X9 A . -)("- 2 )/ 2 

v ; A 2 + 1 + (A 2 - l)cosd(a,w) 

where d n is a uniform constant, d is the geodesic distance on (S n ,c), a G S n , and 
A > 0. For such a function, we will write u a> \ = up- 
Let 

Then it is well-known that the critical points of Jk on S uniquely correspond with 
the solutions of the problem of prescribing scalar curvature function K on S n . As 
shown by A.Bahri (see [15]) and O.Rey (see the appendix D of [17]) that every 
nontrivial u F is non-degenerate on the orthogonal part E 1 - of the subspace E of 
H\S n ) with 

E = span{D 1 , D n , D x }. 

Here Di is the Lie derivative of u F along the conformal vector field generated by the 
orthonormal complement of a G S n (so {a, ei, ...,e n } forms a basis of R n+1 ) and 
L> A is the Lie derivative of up along the scale A. 

Now, we discuss the bifurcation at the point up if H = t.h with t being a small 
scalar and h being a suitable function. We assume without loss of generality that 
the C 2 norm of h is not bigger than one. Because, for t = 0, uf is a non-degenerate 
critical point of on E L , we can use the implicit function theorem to find a 
positive constant e such that for |t| < e there is a solution curve u(t) = Up(t) of 
dJK\exp(E ± )( u ) = f° r an y smooth function ft,. Here exp is the exponential map of 
i/ 1 (S' n ). Now we will find the conditions on h such that 

dJ K \ex P (E){u(t)) = 0. (£>*) 

This means that -u(t) is the solution of the problem of prescribing scalar curvature 
function K = n(n— 1) +t.h. The easy condition for solving (D*) is another (higher 
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order) non-degenerate condition for using the implicit function theorem again. Let 

j(t,y) = J K (t)(u(t,y)). 

Here y = (a, A) G S n x R + are the parameters for u F and we write u(t,y) = u F (t). 
Note u(0,y) = u F (0) = u F . Then we are looking for the conditions for h such that 
u(t, .) is a critical point of j(t, .) for each small t. Note, by the Lagrange multiplier 
method we find that there are some constants h, l n +i such that 

dJ K (u(t, y)) = hDi + ... + l n D n + l n+ iD x . 

In this form we are in the early consideration before Theorem 2. Therefore j y (t, y) := 
d y j(t,y) = is equivalent to U = for all i = l,...,n + 1. By construction, 
j S n Up +1 dfx and we have j y (0, y) = 0. Hence we obtain by using Taylor's expansion 
of j at t — that 

jy(t,y)=tdtjy(0,y)+0(t 2 ). 

So, if our functions d t j y (Q, y) has simple zero at y — y , then we get a solution u(t, y) 
for the problem of prescribed scalar curvature K = n(n — 1) + t.h by the standard 
bifurcation method. In general, one should use the Thom-Mather theory to study 
the zero set of j y (t,y). We leave this open at this moment and perhaps we will 
return to it later. Because up(0,y) is a critical point of Jk at t = and u F (t,y) is 
a curve in E we obtain that 

/ u^(u F ) t (0,y)dfx = 0. 
Js n 

Now, by a direct computation, we find that 

dtj y (0,y) =d y d t j(0,y) 

= f sn hup(u F ) y dfi + d y f sn Up(u F ) t (0, y)dfx, 
= f S n hup(u F ) y dii. 

Hence, we get the following result. 



Theorem 8.2.2 If {u = u F = Ua,x\ is a family of non- constant Yamabe solutions 
and h is a smooth function on S n such that for the orthonormal basis {a, e±, e n }, 
the vector-valued function 

{ hu p u w dfi; w = A, ei, e n } 

has a simple zero or other standard zero point (a , Ao) in the bifurcation theory 
or Catastrophe theory, then there is a positive constant e such that the problem of 
prescribed scalar curvature function n(n — 1) + t.h has a solution, which is near to 
u a ,x f or eac h small t : \t\ < e. 
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We remark here that one can also introduce higher order non-degeneracy if 
J s „ hu N u y dfi is trivial. And the other remark is that this kind of argument can 
also be used to study bifurcation in Nirenberg's problem. The result is 

Theorem 8.2.3 Suppose up = \log\Jp\ is a nontrivial, where F is a conformal 
transformation S 2 . If h is a smooth function such that J s2 he 2uF (uF) y dfi has a 
simple zero, where y is the parameter for the conformal group on S 2 , then there 
is a small positive constant e such that Nirenberg's problem of prescribed gaussian 
curvature 1 + t.h has a solution for each small t : \t\ < e. 

For the detail of our proof of this Theorem, one may see [53] . For more about 
the Yamabe problem, one may see Th. Aubin's book [3] . 
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